Submitted to the Astrophysical Journal 

Preprint typeset using I^IgX style emulateapj v. 10/09/06 



ELECTRODYNAMICS OF MAGNETARS IV: 
SELF-CONSISTENT MODEL OF THE INNER ACCELERATOR 
WITH IMPLICATIONS FOR PULSED RADIO EMISSION 

Christopher Thompson 

CITA, 60 St. George St., Toronto, ON M5S 3H8, Canada 
Submitted to the Astrophysical Journal 

ABSTRACT 

We consider the voltage structure in the open-field circuit and outer magnetosphere of a magnetar. The 
standard polar-cap model for radio pulsars is modified significantly when the polar magnetic field exceeds 
1.8 x 10 14 G. Pairs are created by accelerated particles via resonant scattering of thermal X-rays, followed by 
| the nearly instantaneous conversion of the scattered photon to a pair. A surface gap is then efficiently screened 

by creation, which regulates the voltage in the inner part of the circuit to 10 9 V. We also examine the 
electrostatic gap structure that can form when the magnetic field is somewhat weaker, and deduce a voltage 10- 
30 times larger over a range of surface temperatures. We examine carefully how the flow of charge back to the 
star above the gap depends on the magnitude of the current that is extracted from the surface of the star, on the 
' curvature of the magnetic field lines, and on resonant drag. The rates of different channels of pair creation are 

determined self-consistently, including the non-resonant scattering of X-rays, and collisions between gamma 
rays and X-rays. We find that the electrostatic gap solution has too small a voltage to sustain the observed 
pulsed radio output of magnetars unless i) the magnetic axis is nearly aligned with the rotation axis and the 
light of sight; or ii) the gap is present on the closed as well as the open magnetic field lines. Several properties 
of the radio magnetars - their rapid variability, broad pulses, and unusually hard radio spectra - are consistent 
with a third possibility, that the current in the outer magnetosphere is strongly variable, and a very high rate of 
+3 , pair creation is sustained by a turbulent cascade. 

Subject headings: plasmas — radiation mechanisms: non-thermal — stars: magnetic fields, neutron 
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1. INTRODUCTION 



Magnetars are persistent sources of broad band, non-thermal radiation (see Woods & Thompson 2006 for a recent review). 
The discovery of pulsed radio emission from XTE J1810-197 and IE 1547.0-5408 (Halpern et al. 2005; Camilo et al. 2006, 
2007a,b) adds to the body of evidence identifying them as isolated, non-accreting neutron stars. This radio emission is peculiar 
, in several respects: its brightness fluctuates strongly on a timescale of minutes to days, and peaks above a frequency of 100 GHz. 
The detection of pulsed radio emission from two magnetars, some 15 percent of the presently known population, appears to be 
a consequence of the surprisingly large pulse duty cycle. The known magnetars spin with typical periods of several seconds, 
and their open field lines form a narrow bundle close to the star. For this reason, it had generally been expected that their radio 
emission would be strongly beamed. 

The output of a magnetar in the hard X-ray band can greatly exceed the spindown power of the star (Kuiper et al. 2004, 2006), 
and so most of the underlying dissipation is probably concentrated in the closed magnetosphere (Thompson & Beloborodov 2005; 
Baring & Harding 2007). The same energetic argument does not apply to the radio output, which is smaller than the spindown 
r> | luminosity. The radio emission zone could, therefore, be restricted to the open magnetic field lines. 

It has also been suggested that relatively strong currents flow within the outer magnetosphere of a magnetar, being sustained 
by the outward transfer of twist from an inner zone where it is injected by starquakes. This redistribution of the twist causes the 
poloidal field lines to expand slightly, thereby increasing the open-field voltage (Thompson et al. 2002). The effect depends on 
the extraction of a tiny fraction the helicity that is stored in the inner magnetosphere and stellar interior. The radio brightness of 
the magnetar XTE J 1 8 1 0- 1 97 is observed to correlate strongly with the spindown torque (Camilo et al. 2007a), a behavior which 
is consistent with this type of structural change. Indeed, the torque variations measured in XTE J18 10-197 are comparable in 
magnitude to those seen in the Soft Gamma Repeaters, which are much more variable in the X-ray band (Woods et al. 2002). 

The injection of twist into the outer magnetosphere inevitably triggers current-driven instabilities and strong fluctuations in 
the current. A detailed model of pair creation in a dynamic open magnetic flux tube is given in a companion paper (Thompson 
2008; see also Lyutikov & Thompson 2008). The torsional oscillations of the magnetic field cascade to high frequencies if their 
amplitude exceeds a critical value, resulting in a very high rate of pair creation. 

In this paper, focus on the case where the current that is drawn from the magnetar surface varies only slowly on the timescale 
Rns/c, where R^s is the stellar radius. Our main goal is to work out the structure of a pair-creating diode forming at the 
surface when the magnetic field is stronger than the QED value Bq = 4.4 x 10 13 G. The voltage across the diode is influenced by 
several competing effects: the mechanism(s) of pair creation; the charge flow in the outer magnetosphere; and the gradient in the 
corotation charge density pa at the stellar surface. Our solution applies, in principle, to both the open and the extended closed 
magnetic field lines, because it depends only on the flux of charges returning from the outer magnetosphere. 

A detailed review of the basic QED processes leading to pair creation in super-QED magnetic fields is given in Usov & Melrose 
(1996), and a companion paper (Thompson 2008). The dominant mechanism of pair creation involves the creation of gamma 
rays by resonant scattering of thermal X-rays (Kardashev et al. 1984, Daugherty & Harding 1989, Sturner 1995, and Hibschman 
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& Arons 2001a,b). The voltage across a gap can also be regulated by collisions between gamma rays and thermal X-rays (e.g. 
Zhang & Qiao 1998) and, as we show in this paper, by the direct conversion of non-resonantly scattered X-rays to pairs. 

One of our main results is that a diode cannot be self-consistently be maintained at the surface when the magnetic field 
#ns > 4Bq = 1 .8 x 10 14 G. In this case, the scattered photon is typically above threshold for pair creation, and there is essentially 
no delay between the creation of a gamma ray and its conversion to a pair (Beloborodov & Thompson 2007). We describe a 
circuit solution that has an inward e discharge, within which the voltage that is screened to a low value 10 9 V. The stability 
of this solution is also investigated, and is shown explicitly to depend on a delay between gamma-ray emission and pair creation. 

A surface gap is possible when Z?ns < 4Bq and the current density is larger than pac. We construct a detailed model including 
all the relevant modes of pair creation. We focus on space-charge limited flows, and explain why charges of both signs can be 
supplied continuously from the surface of a magnetar. Even through heavy atoms are tightly bound into molecular chains (Medin 
& Lai 2007, and references therein), ions can be liberated from the magnetar surface by photoionization and knockout processes. 

The equilibrium gap voltage is calculated as a function of the X-ray black body temperature 7bb of the stellar surface. It is 
gradually reduced as A:B7t,b rises above ~ 0.15-0.2 keV, and non-resonant scattering and photon collisions supply increasing 
amounts of pair creation within the gap. Pairs that are created just outside the gap continue to lose energy by resonant scattering, 
and the resulting increase in the pair multiplicity is derived. We explain in detail how the polarization of a pair-creating cloud 
that is formed just outside the diode will create enough returning charges to compensate the gradient in corotation charge density 
outside the gap. 

An important subtlety involves the effect of cyclotron drag at ~ 10 stellar radii; we show that outside this zone the correct space 
charge is obtained if the charges of sign opposite to pa flow outward sub-relativistically. This solution breaks down beyond a 
critical angle from the magnetic axis, where charges of both signs can be decelerated to trans-relativistic speeds. As a result, we 
point out that there is a selection effect against observing pulsed radio emission from X-ray bright neutron stars with spin periods 
shorter than ~ 0.3 seconds. The reflection of gamma rays back toward the star via resonant scattering can have a similar effect at 
much longer spin periods. 

The plan of the paper is as follows. The origin of the open-field voltage is reviewed in § [2j and the two basic types of circuit 
solution - with and without a diode - are described. The voltage of a surface diode in a magnetic field < 4Bq is derived in $3] 
and the effects of cyclotron drag are investigated in §|4] In section[5]we explore in detail the plasma flow on the open field lines 
when Z?ns > 4Z?q and a diode is absent at the surface of the star. The paper closes with a summary of how the different circuit 
solutions described in the paper are constrained by the luminosities of the radio magnetars. The Appendix compares the effects 
of field-line bending and general-relativistic frame dragging on the gradient in corotation charge density. 

2. VOLTAGE GENERATION 

In the absence of pair creation, a very large voltage develops on the open magnetic field lines of a rotating neutron star. When 
the effects of screening by pair creation are included, the magnitude and time dependence of the voltage are believed to depend on 
several factors (e.g. Sturrock 1971; Ruderman & Sutherland 1975; Scharlemann, Arons, & Fawley 1978). A vaccum gap could 
form with a relatively large and variable voltage when the corotation charge density is positive, if the ions are tightly bound in 
molecular chains to the neutron star surface. Although magnetars are indeed expected to have condensed surfaces in the absence 
of a hydrogen or helium layer (Medin & Lai 2007), it is p laus ible that photodissociation and spallation process will liberate ions 
at a sufficient rate to supply the spindown current (Section l2~3"T l. We therefore assume a continuous charge flow from the magnetar 
surface, for either sign of pcj- 

The voltage that develops in this case has generally been though to be regulated by the gradient in pa away from the neutron 
star surface. The dominant contributions to this gradient come from general relativistic frame dragging (Muslimov & Tsygan 
1992) and magnetic field line curvature (Scharlemann et al. 1978). The second effect may well be dominant in magnetars, given 
the strong evidence for a non-dipolar field structure in the light curves of SGR flares (Thompson & Duncan 2001). 

The main assumption here is that the current that is drawn from the neutron star surface can be supplied very nearly by a flow 
of outgoing and returning charges, 

J — P>v> + P<v< — (p>-p<)c (1) 

In our basic one-dimensional model, the current density J is nearly aligned with the background magnetic field, and p^ and p^ 
are the densities of outgoing and ingoing charges. The flux of ingoing charges is provided by pair creation within the circuit. 

We consider two basic circuit solutions. In the first case, the pair creation is distributed smoothly throughout the circuit, so as 
to maintain the local charge balance 

Pt> (z) + p< (z) = pc] (z) . (2) 

Here pa is the corotation charge density. The second basic type of circuit solution involves a plasma-filled gap of height h, which 
sits at the base of the circuit. Charges moving outward through this gap gain enough energy to emit gamma rays that are able to 
convert to pairs at z > ft, thereby creating a dense polarizable plasma. The gap voltage can also be limited by pair creation by 
ingoing charges, and by modes of pair creation that are distributed throughout the volume of the gap 
In a Newtonian gravitational field, the corotation charge density would be given by (Goldreich & Julian 1969) 

^ = ~V.(txB) = -^. (3) 

4irc lire 



This result is modified by relativistic effects: in particular, the rotation of the neutron star drags the frame of a zero angular- 
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momentum observer with an angular velocity us = kSI, where the coefficient is given by 1 

n j «NS="^ T =0.15/ NS .45^ N S6- ( 4 ) 

Here 7ns is the moment of inertia of the neutron star, and 7?ns its radius. This has the effect of reducing the corotation charge 
density, 2 

Pa = Pgj[1 -«(>•)]• (5) 

The normalized charge density pa/B increases with radius along a straight magnetic field line. 

If the higher order multipole components of the magnetic field are comparable in strength to the dipole near the neutron star 
surface, then a stronger voltage can be produced than purely by the frame dragging effect (Barnard & Arons 1982; Asseo & 
Khechinashvili 2002). As we show in Appendix lAl the leading contribution to the field line curvature comes from the octopole 
component of the magnetic field, and the quadrupole component to second order, when these two components are treated as a 
perturbation to the dipole, 

d(p G3 /B) p G3 . , ~ PG3 ( r \~ 3 

B — = sinxcos0 c = K 3 . (6) 

dl Rc Rns \Rns J 

Here Rc is the curvature radius and the coefficient kj, is expressed in terms of the surface amplitudes £2.3 of the quadrupole and 
octopole via 

£3= (^3 + e 2 ^^J sinxcos0 c - (7) 

In addition, \ is the angle between the magnetic moment fj, and the rotation axis fi, and <f>c the angle between the plane of 
curvature of the magnetic field lines and the /x-fi plane. 

A subtlety associated with this model, which has not been fully explored, is that the value of J that is demanded by the outer 
magnetosphere need not closely match an outward charge flow of density pcj and speed c. When the current density is smaller in 
magnitude than pac and has the same sign, the open-field voltage remains modest close to the star, 3 and the mean flow speed of 
the outgoing charges is smaller than the speed of light (Mestel et al. 1995; Shibata 1997; Beloborodov 2007). We therefore focus 
here on the case where J is larger in magnitude than p G sc, allowing for the possibility that J may have the opposite sign. Such an 
imbalance could naturally result from a multipolar structure in the external magnetic field of the magnetar; and from more subtle 
effects such as resistive instabilities in the current carrying flux tube or more exotic effects such as the loss of magnetic helicity 
from the closed magnetosphere onto the open field lines. 

2.1. Voltage Solution I without an Inner Gap and Pair Multiplicity ~ 1 

We now describe a simple model for the inner part the open field circuit, within which pair are created at the minimum rate 
needed to balance the gradient in p G3 along the magnetic field. This zone extends from the neutron star surface out to a maximum 
radius which is determined self-consistently by the stability of the solution ( jj5.ll ). The range of radii considered here is 
much broader than the width of the open-field bundle, i.e. -Rns ^ ^NS^openC^NsX where 6* op en('") = (S!r/c) 1//2 . The circuit can, 
therefore, be approximated as a thin tube. 

We take the flow to be steady, with all particles moving relativistically and charges of opposite signs moving in opposite 
directions. Implicit in this setup is the assumption that one particle in each newly created <? ± pair reverses direction nearly 
instantaneously. The flows of ingoing and outgoing charges satisfy the usual continuity equations, 

d 

c — (n > S±) = n±S±, (8) 
ar 

-c^-(n < S±) = h±S±, (9) 
dr 

where ri± = (n e - + h e +)/2 is the rate of creation of electrons and positrons per unit volume, and S± (eq. IIA16I ) is the cross sectional 
area of the open field-line bundle at radius r. The gradient in the total charge density p = Pt> + P< is obtained from the sum of 
equations © and ©. One finds 

d 

c — (pSi_) = 2en±S±sgn(p G] ), (10) 
ar 

where e is the magnitude of the electron charge. 

As is shown below, the electric field Eh = E • B is greatly reduced compared with the unscreened value in the absence of pair 
creation. The small Eu is sustained by a small fractional departure of p from the local corotation charge density. Therefore, 

p(r)~p G1 (r). (11) 

1 Throughout this paper we use the shorthand X = X„ X 10", where the value of quantity X is measured in c.g.s. units. 

2 The quantity on the left-hand side is, strictly, the product of the lapse function a(r) and the proper charge density, but this subtlety can be neglected in our 
present treatment because B and the current density apQjv both satisfy the same conservation equation. 

3 This statement does not apply to the outer magnetosphere, where part of the open-field bundle crosses the surface where ■ B vanishes; in other words, an 
outer gap can still be present even in the absence of an inner gap. 
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Associated with this charge density is a net outward flow of particles moving at v : 
The continuity equation ( TT0l > can then be written as 



c, 



-\PQ}\S±C = 

e e 



I dNg 
2 dr 



(12) 



(13) 



A specific example in which pair creation is driven by resonant scattering is described in $5] 

In this toy model, the net current along the open field-line bundle is in a steady state, 

1^+1^ = 1 = const, (14) 
where 1^ = cp t> S± and /<, = -cp < S±. The strength of the current is determined by conditions in the outer magnetosphere, and 
does not influence the qualitative nature of the solution. In particular, current-driven instabilities in the outer magnetosphere 
are assumed to occur on timescales much longer than the light-crossing time of the inner circuit, so that changes in / occur 
adiabatically. (This assumption will break down if the outer magnetosphere is strongly dynamic, so that a turbulent spectrum of 
modes is established by non-linear couplings between the waves: Lyutikov & Thompson 2005; Thompson 2008.) 

Using p ~ pcj we also have 

/>-/<=/gj- (15) 

Equations (TFfl i and ( fTBI l imply 

t 1 (j t \ dI< 1 dIa ntz\ 
I<=~(I-Igj), —r- =—x—r~ ■ (1") 



2 dr 2 dr 

The outward current density can therefore differ significantly from /qj if / does, 

I> = ^(/+/gj). 



(17) 



We allow for the possibility that a gap may form beyond some outer radius This radius is determined, in the case of a 
magnetar, by the condition that resonantly scattered X-ray be able to convert directly to electron-positron pairs, and corresponds 
to the surface where B = 4Bq (©. Beyond this gap, a dense polarizable plasma is present and the solution continues to hold in a 
time-averaged sense. 

2.2. Voltage Solution II with an Inner Gap and Pair Multiplicity 3> 1 

In a super-QED magnetic field, the conversion of gamma rays to pairs generally occurs over a small distance compared with 
the width of the open magnetic field lines. We therefore focus on a gap with height smaller than ~ # pen^NS- Particles flow 
relativistically through nearly the entire gap width, and so p< and are taken to be constant. The electrostatic potential satisfies 
the equation 

<9 2 $ 

—— = -4 7: (p t> + p < -p C] )_ (18) 
Close to the surface of the neutron star, the corotation charge density has an approximately linear dependence on height, 



PGi = Pgs{Rns) 1 +A 



(19) 



The dimensionless quantity A can be expressed in terms of the frame-dragging coefficient k (eq. J4)) and the octopole-quadrupole 
coefficient K3 (eq. Q) via 

A ~ 3k(R ns ) + K3. (20) 

(see Appendix lAl for further details). The pair creation rate is assumed to be zero at at height z < h above the magnetar surface; 
and to become large at z > h. In other words, the density of pairs that are created outside the gap is assumed to become much 
larger than the change in pes je across the gap, within a small distance z-h<sch from the upper boundary of the gap. 

Integrating eq. (fT8l with respect to z with the boundary conditions $ = t/$ fdz = at z = and t/$ jdz = at z = h implies a 
unique value of the charge density inside the gap, 



P = P> + P< = PCS z 



2 



(21) 



The gap solution is 



P> = ~ 



P< = =■ 



J h 
c +PG3 {2 

J (h 
- C +PCS {2 



d$ h 2 
E z (z) = -— = 2TrAp GS (R NS )— 
dz «NS 



*(z) = 3» g 



gap : 



7T h 3 
-W A PG1( R NS) — - 



(22) 
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These expressions for E z and $ are very similar to the diode solutions constructed by Arons & Scharlemann (1979) for a space- 
charge limited flow in a planar geometry, the main difference being that a returning flow of charges is included. 

The establishment of the screened state at z > h typically requires a very modest additional potential drop A<£> compared with 
the gap voltage d22b . Outgoing charges of the sign opposite to pa must reverse direction so as to maintain p< + p> = pair)- 
The outgoing charge density has both positive and negative components, pi> = p^ + p^,, with \p^\, \p^ \ 3> |p<|. The density of 
ingoing charges is just sufficient to compensate the difference between J and pac, 



p<(z) = -~ 



J 

--PGi(z) 

c 



When the pair creation rate is high at z > h, there is a small step in the ingoing charge density aXz = h, 



p < (h+e)-p4h-e) = ~ 



PGiW-pci ( - 



A f h 

4 l%s" 1 PGJ(%S) - 



(23) 



(24) 



We work out a specific example of such a gap in 5(3] where pair creation is dominated by the resonant scattering of thermal 
X-rays. When the surface magnetic field is ~ 10 14 G, particles crossing the gap acquire an energy e|$ gap | ~ 10 3 - 10 4 m e c 2 
(depending on the surface temperature). 

2.2.1. Sustaining the Low-Voltage Circuit Solution Outside the Gap 

The gap solution ( 122b depends on the ability of the circuit to maintain a configuration with a relatively high particle density 
outside the gap, 

>Gj| 



and an approximately charge-balanced state with a low voltage, 



P>(r) + p<i(r) ^ Pair)- 



(25) 



(26) 



In the inner part of the circuit, this second condition can be satisfied if some fraction of the 'negative' charges (with a sign 
opposite to pcj) reverse direction. This is possible because the outgoing primary charges that flow through the gap will create a 
dense cloud of secondary pairs as they lose energy to resonant drag. The energy spectrum of these secondary particles extends 
far below the gap voltage. As a result, only a small additional voltage is needed to create the required number of ingoing charges. 
Indeed, the photons that are resonantly scattered by the outgoing primary charges will continue to create pairs to a much l arge r 
distance from the star, where B~4x 10 12 G and the energy of the secondary pairs is ~ 10- 100 times smaller than e$ gap ( j34.lt . 

Screening of the electric field near the outer gap boundary is straightforward in the solution (l22l . because £ji — > inside the 
upper boundary of the gap even in the absence of pair creation. The minimal pair density that is needed for screening is higher 
when Eii at the first point of pair creation remains comparable to the peak value in the circuit (e.g. Shibata, Miyazaki, & Takahara 
1998, 2002). This is the case in the vacuum gap model of Ruderman & Sutherland (1975). 

At greater distances from the star, where the cyclotron energy of an electron is ~ ^B^bb, mildly relativistic electrons and 
positrons feel a very strong drag force (e.g. Thompson & Beloborodov 2005). In this zone, an outgoing charge can reverse 
direction only in the presence of a relatively strong En, For the moment, we will neglect the effects of radiation transfer and 
assume that a large fraction of the outflowing X-rays are unscattered. Near the magnetic pole, an electrostatic potential 



ed$ = - 



- eJ B|| fl fr = 8.0xl0 2 J B- 1 ^ 5 <,6 



(dL x /duJx)duJx 
10 35 



ergs s 



VkeV 



-2/3 



MeV 



(27) 



just cancels the outward force of resonant scattering on the electrons. The dependence of <f> on radius is plotted in Fig. Q] assuming 
a blackbody X-ray spectrum. When the magnitude of this potential is much larger than the kinetic energy of the secondary pairs, 
it is not possible for the circuit to maintain a steady state in which only a fraction of the secondary pairs reverse direction. 

A second steady-state circuit solution can, however, be constructed at a somewhat smaller voltage. Here we drop the assumption 
that all the charges move close to the speed of light. An electric potential just slightly smaller in magnitude than eq. $T7\ will 
allow the negative charges to flow outward at a speed somewhat below the speed of light, even while the outgoing positive 
charges remain highly relativistic. The proportions of the current J+, J- that are carried by the positive and negative charges can 
be left arbitrary, and related to the respective space charge densities via J+ ~ p+c, J- = p-V— Demanding that J+ + J-=J and that 
condition ( f26l > be satisfied, we find 

J~Pgjc\~ \J-\ 



-= 1 + 



\J- 



J+ - pac 



(28) 



This drift speed is indeed smaller than the speed of light, since we are assuming that the total current density J > pac (otherwise, 
the negative charges are not required to maintain local charge balance within the circuit). The transition between the two circuit 



solutions occurs at the radius where eEur 



In 



e c , where 7 m i n is the minimum Lorentz factor of the secondary pairs. 



The fluctuations in the current density are not small when the twist on the open magnetic field lines is dynamic, and there is 
a significant coupling between the ingoing and reflected torsional waves. The resulting cascade generates very strong current 
fluctuations, and most of the energy in the fluctuating toroidal magnetic field is deposited in the particles (Thompson & Blaes 
1998). In this case, the creation of pairs does not depend on large scale electric fields, and very high pair multiplicities are 
possible. This type of circuit solution is examined in Thompson (2008). 
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FIG. 1. — Electric potential gradient (given by eq. 1271 ) that is required to cancel the radiative force exerted by thermal X-rays on a stationary electron or 
positron. A blackbody spectral distribution is assumed, with emission from the surface of the neutron star (Rns = 10 km) in a single polarization mode. The 
different curves correspond to X-ray luminosities (1,3, 10,30, 100) X 10 33 ergs s (from bottom to top). If a surface gap emits dense clouds of outgoing pairs, 
then an electric potential close to this value will develop beyond the radius ~ (10— 15)i?Ns where the potential drop e(r/3)d&/dr is larger than the minimum 
energy of the outgoing charges. This relatively weak electric force allows the outflowing plasma to maintain the local charge balance p = pgj('") in this zone. 



Finally, we note that our gap solution makes sense only if the outgoing and ingoing charges within the gap have opposite signs. 
One sees from eq. (1221 that this is possible only if J is larger in magnitude than pgj(Rns)c When J is opposite in sign to pcjc 
but smaller in magnitude, then the current can be supplied by an inward subluminal drift of the corotation charge, but only if 
there is a source of charge in the outer magnetosphere, e.g, associated with an outer gap, or a diversion of current across the 
magnetospheric boundary. 

2.3. Supply of Positive Charges from the Surface of a Magnetar 

The current on the open magnetic field lines is supplied mainly by electrons when the spin vector and dipolar magnetic field 
are aligned (fl ■ B > and pcj < 0); and by positive charges (ions or positrons) in the case of anti-alignment (fi • B < and 
Pa > 0). It has been suggested that the observational manifestations of magnetars could be substantially different in the two spin 
orientations, because heavy ions would be tightly bound to the surface of the star in long molecular chains (Zhang & Harding 
2000). Although the binding energy of atoms heavier than oxygen is indeed too high to allow thermionic emission (Thompson et 
al. 2000), two mechanisms of creating unbound ions are available. 

First, downward-moving relativistic electrons trigger a cascade of pairs and gamma rays below the condensed surface. Gamma 
rays with energies of ~ 20-30 MeV are created in such a cascade, and are able to knock out protons and neutrons from heavier 
nuclei. On average ~ 0.1(7/10 3 ) proton is created per seed relativistic electron of Lorentz factor 7 (Beloborodov & Thompson 
2007). The pairs created in a surface gap (where B < 4-Bq) have typical Lorentz factors 7 ~ 10 4 (§|3), an d me magnitude of the 
returning charge density can approach pa if there is a significant mis-match between J and pac. In this case, knockout appears 
marginally effective at supplying a Goldreich- Julian flux of particles. However, the energy of the returning charges is regulated 
to 7- 10 3 ifB NS >4B Q (©. A gap may therefore form above the condensed surface with a voltage regulated to the value that 
gives about ~ 1 proton knocked out per returning charge. This voltage (~ 10 10 V) remains too small for the accelerated particles 
to emit curvature photons in the gamma ray band. An additional population of energetic charges could possibly be created in the 
outer magnetosphere, with Lorentz factors high enough that returning to the star they suffer negligible resonant drag. (The net 
kinetic luminosity of these returning charges need only reach ~ I0 4 m e c 2 lcj/e, typically ~ 10~ 3 of the spindown luminosity of a 
magnetar.) 

Second, ions can be photodissociated from the molecular chains in the presence of a strong flux of 10-100 keV photons. 
Most active magnetars are observed to emit such hard X-rays at a rate Nx ~ 10 42 - 10 43 s -1 (Kuiper et al. 2004; Mereghetti 
et al. 2005). The observed spindown of the neutron star implies a much smaller particle flux along its open field lines: a 
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spindown luminosity L sc j corresponds to a total Goldreich- Julian flux Na ~ 1 X lO 31 ^',^ s This flux arises from a fraction 

ttRns/cP = 2 x 10~ 5 (P/6 s) -1 of the neutron star surface. If even one part in a million of the hard X-ray flux is reflected toward a 
part of the stellar surface that contains the polar cap region, then the photodissociation rate in the polar cap could be high enough 
to sustain Ni on ~ Ncs- To put this number in perspective, it will be recalled that a twist angle of ~ Acf> radians in the closed field 
line region of a magnetar supplies an optical depth ~ A<fr at the electron cyclotron resonance. 

Ions themselves are not efficient radiators, so the question then arises as to how pair creation can be spawned in the open field- 
line region when pa > 0. The same issue is encountered in electrodynamic models of radio pulsars (e.g. Cheng & Ruderman 
1977). The simplest answer is that a state in which the 'excess' corotation charge density pa - Pes ( ec l- 0) is supplied by pair 
creation is self-sustaining, just as it is in the twisted magnetosphere of a magnetar (Beloborodov & Thompson 2007). The details 
of how this happens are the subject of §§0and[5] 

3. GAP IN A MAGNETIC FIELD B < 4B Q 

The strength of the surface magnetic field has an important influence on the mechanism of pair creation (Beloborodov & 
Thompson 2007; Thompson 2008). Magnetars are bright thermal X-ray sources, so that relativistic charges produce pairs in 
copious amounts by the resonant scattering of the X-rays. When B > 4Bq, some resonantly scattered photons are created above 
the threshold for pair creation, and there is essentially no delay between the positions of scattering and pair conversion. This 
leads to effective screening of the component of electric field parallel to B. 

The critical field strength for direct pair creation of this type is obtained easily by considering the Landau de-excitation of an 
electron (or positron) and the energy of the resultant gamma ray. In a frame in which the excited particle is at rest, the gamma 
ray energy is maximized when it is emitted orthogonal to B. In this case, the electron receives no recoil, and energy conservation 
gives (1 +2B/Bq) 1 / 2 ot c c 2 = E 1 + m e c 2 . Requiring £ 7 > 2m e c 2 corresponds to B > 4Bq. 

In weaker magnetic fields, direct pair creation by non-resonant scattering is possible at a much reduced cross section. The 
target photon must have an energy larger than the resonant energy hui' le;i = (B/BQ)m e c 2 in the rest frame of the scattering charge. 
We consider the e ffects of this process in 33.3. 1 1 Collisions between gamma rays and thermal X-rays can also be important, and 
are examined in 33.3.21 The resulting zonal structure is summarized in Fig. [2] Many more details and references about pair 
creation in ultrastrong magnetic fields are given in Thompson (2008). 

In this section we consider the case where B < 4Bq at the surface of the star, so that the conversion of resonantly scattered 
gamma rays to pairs occurs only adiabatically. First, the gamma ray propagates in the curving magnetic field and converts to 
bound positronium; and, second, the positronium atom is dissociated when one of the constituent charges resonantly scatters an 
X-ray or absorbs an infrared photon. 

3.1. Resonant Scattering of Thermal X-rays 

The mean free path for the creation of gamma-rays, and then for the dissociation of the positronium atom, is obtained from the 
rate of resonant scattering. An X-ray of frequency 

u' x = 1 {\-pp)uj x = — (29) 
m e c 

in the rest frame of an electron or positron is scattered with an enhanced cross section. Here p is the direction cosine of the target 
photon with respect to the magnetic field (defined so that p > and j3 > when the photon and charge move away from the star). 
The target photon is strongly aberrated into a direction nearly parallel to B when the charge moves relativistically (f3 = v/c ~ 1), 
so that 

o»(a£)~— sL' x -— ) (30) 
m e c \ m e c J 

(Daugherty & Ventura 1978). We focus on the case where the target photons are emitted spherically from the surface of the star, 
with a temperature 7bb and a single polarization state (e.g. Silantve & Iakovlev 1980). The effects of light bending are neglected, 
so that p uniformly fills the interval I > p > /i m ; n , where 

Mmin = (l-4sA 2 ) 1/2 - (3D 

One has p m [ n = at the surface of the star, and /i m ; n ~ R^ s /2r 2 at r 3> R^s- 
A charge moving radially outward with a Lorentz factor 7 > 3> 1 scatters X-rays at the rate (Sturner 1995) 



r res _ a / Q bb n ,, . -.2-2 



2(B/B Q ) 

Here 



]n(l-<f a '>) 



c 



(32) 



„ _ Hlj x B/B q 
w> = 7— 7^ = — n n r, (33) 

KB •'bb 7> (= 'bb( 1 ~ Mmin) 



is the dimensionless frequency of the target photon at p, = p m i n - We sometimes use the more compact notation 

2>= n =2.3 xlO^l-Mmtor 1 Sis f^f^r ) (34) 

(1-/W)7> V°- 5keV 
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charge flow 



12 



B~4x10 G 



= 4B 



QED 



Zone 1 



Zone 1 : 
Direct pair creation by 
resonant scattering 

Zone 2: 

Indirect pair creation by 
resonant scattering 

Zone 3: 
Secondary pairs below 
energy threshold for 
continuing pair creation 

FIG. 2. — Zones of pair creation on the magnetic field lines extending to a large distance from the magnetar surface. The direct conversion of resonantly 
scattered photons to pairs is restricted to the innermost Zone I, where B > 4Bq. Where the magnetic field is weaker, resonantly scattered photons convert to pairs 
only after propagating a macroscopic distance. When a dense cloud of pairs forms in a gap close to the surface, large scale electric fields are screened, but new 
pairs can continue to form as the outflowing particles lose energy to resonant drag, out to the radius where 8~4x 10 12 G. Secondary modes of pair creation are 
important inside the gap when the surface temperature > 0.2 keV. Target photons that are non-resonantly scattered with a reduced cross section can convert 

directly to pairs in Zone 2, and collisions between resonantly scattered photons and X-rays are important at somewhat higher surface temperatures. See the text 
for further details. 



Let us evaluate expression d32l i at the surface of the star. An outoing charge experiences a large number of scatterings if the target 
photons are drawn from the black body peak (u}[> ~ 3), 



^T™(R NS ) = 3.9xl0 3 B^ 15 
c 



^0.5 keV 

An analogous expression is obtained for the scattering rate of an ingoing charge 



ln(l-e"^) 



pres 
1 < ' 



O/Obb ( B 



Br 



In 



1 - exp ( 



1 — exp 



( B/Bq \ 

V 2 7< ,e bb y 



c 



which reduces to 



B/B 



Q 



far from the stellar surface. Here = w^(e w<1 - 1) , 



w<(7<)= — , 

7< 



(35) 



(36) 



(37) 



(38) 



ELECTRODYNAMICS OF MAGNETARS: INNER ACCELERATOR 



9 



Surface Gap (B NS < 4B Q ) 



r = R 
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FIG. 3. — Structure of a gap at the surface of a neutron star with a polar magnetic field Bns < 4Bq. Some key modes of pair creation are labelled. Outgoing 
particles resonantly upscatter thermal X-rays that are emitted from the surface of the neutron star. When the surface temperature is lower than ~ 0.2 keV, the 
width of the gap is is set by the conversion of these gamma rays to free pairs off the magnetic field. (The gamma ray first converts to bound positronium, which 
is then rapidly dissociated when one of the charges resonantly scatters an ambient X-ray.) A fraction of the gamma rays within the gap can convert to pairs by 
colliding with an X-ray before reaching the threshold energy for single-photon pair creation, and direct pair creation by non-resonant scattering is also important. 
The gap voltage is lowered by these effects when ^B^bb ;J 0.2 keV. 



and 



0.5 keV / 



(39) 



3.2. Basic Gap Solution 

We now determine the basic properties of the gap solution outlined in 32.21 In this solution, the total current density is allowed 
to exceed pac, with a returning flux of charges of the sign opposite to pes supplying a significant fraction of J. The gap voltage 
is then proportional to the gradient in pa along the open magnetic field lines, 



I* 



gap | 



k dp C ] 3 _ 7r pcj/r 
3 dr 3 /?ns 



(40) 



just as it is in models of pulsar circuits with a steady, space-charge limited flow (Arons & Scharlemann 1979). In contrast 
with models of pulsar polar caps with realistic pair creation (e.g. Hibschman & Arons 2001a), the thickness h of the gap can 
self-consistently be assumed small compared with the width of the open magnetic field lines, h <C #open^NS- 

We start with the simplest solution, in which <£> gap and h are determined by i) the gradient in pes away from the neutron star 
surface; and ii) the degree of curvature of the magnetic field lines (Fig. |3j. This solution is valid if pair creation occurs mainly 
by the adiabatic conversion of gamma rays off the magnetic field, and the charges streaming through the gap feel a negligible 
drag force due to resonant scattering. It turns out that both of these conditions are satisfied if the surface temperature is below a 
critical value, k^T^ < 0.15-0.2 keV. 

The thickness of the gap is determined by summing three quantities, and then minimizing this sum: i) the distance Az for an 
outgoing particle to be accelerated to the Lorentz factor 7 at which it resonantly scatters a background photons; ii) the distance 
Az' for the resonantly scattered photon to convert to a pair; and iii) the distance Az" for one of the created particles to reverse 
direction in the background potential, thereby polarizing the plasma within the gap. This procedure is self-consistent when the 
mean number of scattered photons per outgoing charge satisfies, 



A/?< Ah 



Pa 4/? N s 



(41) 



Here Ap < = p^Qi + e) — p^Qi — e) is the shift in the ingoing charge density across the outer boundary of the gap, and is given 
by eq. (124-b . The right-hand side of eq. d4lT i is small, and so this condition is consistent with weak drag. In this case, the gap 
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FIG. 4. — Possible structure of a gap whose inner boundary is shifted just outside the surface Bns = 4Bq, in the case of a magnetar with a polar magnetic field 
Bns > 4Bq. The existence of such a gap depends on the non-linear development of the two-stream instability within the circuit. 



thickness depends only weakly on the density of target photons, and we can take the energy of the resonantly scattered photon to 
be the maximum possible energy. In a super-QED magnetic field, this is a sizeable fraction of the kinetic energy of the scattering 
charge, 

/•.. (max) = 



1+2B/B Q 



(42) 



In the regime of weak drag, multiple scatterings are rare and a scattering particle attains its kinetic energy jm e c 2 at a potential 
drop A$ ~ 7«i e c 2 /e from the gap boundary. After conversion of the gamma-ray to a pair, one of the two outgoing charges will 



reverse direction after the particles pass through an additional potential drop A$" 
energy and position within the gap is 



7 



7>r 



l~fm e c 2 /e. The relation between particle 



(43) 



where 



In 



A /^ns\ eB NS R NS f _h_\ 3 = 



6 V 



m e c 2 \R N s J 



cap 



h 

^NS 



(44) 



(see eq. 1221 ). The position of scattering is a distance Az from one boundary of the gap, and the position of pair conversion is a 
distance Az" from the other boundary, where 



A£ 

h 



7 



In 



Az" 



7/2 

Tmax 



(45) 



As the thickness of the gap shrinks, so does the voltage and the particle Lorentz factor. The delay between gamma ray emission 
and pair conversion therefore plays an essential role in stabilizing the gap. The conversion distance Az' for the gamma ray scales 
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inversely with 7, and so increases with decreasing h, 

2m P c 2 



Az' 



-Rc 



2 ( 2B/Bq V 



Rc 



(46) 



£ 7 (max)~^ j\l+2B/B Q J 

Here Rc is the curvature radius of the magnetic field lines, which in the case of strong field curvtaure is given by eq. dA9t , The 
equilibrium gap width can be obtained by setting 



f^ = !-(Az+Az' + Az")=0. 

Differentiating eqs. (|45l l and (|46] l with respect to 7, and substituting into eq. ((47) gives 

1 1 _ Az'/h 

(f>'(Az/h) + 24>'{Az"/h) ~ <t>{Az/h)' 
Combining this with the relations cj)(Az/h) = 2(f>(Az" /h) and Az' = h — Az- Az" „ we obtain the solution 



Az~0.38/j; Az"~0.25/!, 



and 



= 17.5 



2B/B Q 
1+2B/B C 



Rc 
h ' 



(47) 



(48) 



(49) 



(50) 



Substituting again for 7 max gives the equilibrium gap width, 



h 
Rns 



2B/B Q 
17.5 \1+2B/B Q 



^gap 



? l/2 
l NS,6 



2.7 x 10 



' A D NS,14 



-1/4 
P 

6s 



1/4 



and voltage 



where 



m e c 



6-5 xlO 3 1/4 1/4 I/2 

J^IHB) B NS,14 A NS,6 



p 

6s 



1/4 



-1/4 



cm, 



PnS / 



3/4 



2B/B Q 
\+2B/B Q 



(51) 



(52) 



(53) 



We now consider the number of photons that are resonantly upscattered by the outgoing charges. The scattering rate is, in fact, 
somewhat higher for the ingoing charges (which see photons moving toward them nearly head on). The maximum energy of the 
resonantly scattered photons is, however, the same for particles moving in both directions, and the distribution of kinetic energies 
within the gap is symmetric between the two directions. 4 The gap width, as determined by the minimization procedure j471 l, is 
therefore independent of the direction in which the scattering charges are moving. 

A simple estimate of the energy of the gamma ray is obtained by viewing the scattering process as resonant absorption followed 
by de-excitation. The absorbing charge decelerates from a Lorentz factor 7 to 71 < 7, which is found from the equation of energy 
conservation, huix + , fin e c 2 = j\Eb, where Eg = (1 + 2B / Bq) 1 > 2 m e c 2 is the energy of the first Landau level. This gives (neglecting 
the small energy Tilox) 



7i 



7OT e C 



: 7 



1 +2 A 



-1/2 



(54) 



The final Lorentz factor of the electron, when it returns to the ground state, equals j t only if the de-excitation photon is emitted 
perpendicular to B in the electron frame. Otherwise there is a recoil effect whose sign depends on the photon direction. If we 
adopt 71 as an estimate of the mean final Lorentz factor, the typical loss in scattering is 5 

A7~7-7! =/ rec0 ii7) 

where 

B 

/recoil = 1 - I 1 +2- 



B 



-1/2 



(55) 



(56) 



The number of scatterings per outgoing particle is distributed over gamma ray energy E- 1 = f KCO i\ini e c 2 according to 



dM. 



dE~, 



d~, 



^fm e cY 



res 



eE, 



(57) 



4 Each outflowing particle scatters X-rays at a modest rate. We therefore neglect the effect of the recoil on the motion of a typical outgoing particle. 

5 It equals 1 7 for B = 4Bq and ~ 0.85 7 for B = 10 15 G. The lost energy is carried by the scattered photon. It is shared between two new particles when the 
photon converts to a pair. Each particle in the created pair therefore has a typical kinetic energy ~ 0.3-0.4 7 m e <p- in the direction of the scattering charge. 
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FIG. 5. — Ratio of drag force due to resonant scattering, and the accelerating electric force, within the gap solution given by eq. (22). The forces on the 
outgoing and ingoing charges are evaluated seperately, as a function of the distance from the injection radius (z = 0. h respectively). The results are displayed for 
surface magnetic fields B^s = Bq and 3Bq. Each set of five curves corresponds to blackbody temperatures kgT^ = 0. 1 , 0.2, 0.3, 0.4, and 0.5 keV (thicker curves 
representing higher temperatures). We choose parameters A = 1, Z?ns = 10 km, Rc = Rns and P = 6 s. 



where we take /j, m i n = at the surface of the star. Substituting eqs. d22"l >. d32i >. d5T1 >. and d52l > gives 



dM- 



dE~, 



= 1.46 



a f Q 



g>|ln(l-e-^>)| 
(z/h)-(z/h)i 



AV2RV2 lo.5keVy 



p 

6s 



1/2 



Rns 



-1/2 



5>|ln(l-e^">)| 
(z/h)-(z/h) 2 ' 



(58) 
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FIG. 6. — Critical blackbody X-ray temperature above which charges flowing through the middle of the gap (22) feel a drag force equal to eE^ , as a function 
of magnetic field in the gap. We choose parameters A = 1, i?NS = 10 km, Rc = Rns and P = 6 s. 



where the normalized energy of the target photon is 



:0.36 



NS.14 



A K NS,6 



0.5 keV J 



P 

6s 



1/4 



-3/4 



7©bb V 7 

The analogous expression for ingoing charges is obtained from eq. ( 1371 ) for the scattering rate r^; one finds 



dM- 



dE~, 



ln(l-e" 



>/2) 



ln(l- 



>) 



1 



aLA/f. 



(59) 



(60) 



where we have made use of the relations ui^ = icD^ for /j, m \ n = (eqs. 1331 , l38ll ). 

A useful measure of resonant scattering within a gap is the relative magnitude of the electric force and the drag force due to 
scattering, 



drag 

eEn 



= /r 



dM~ 



recoil 



d\nE 



(61) 



This quantity is plotted in Fig. [5]as a function of blackbody temperature for polar magnetic fields Z?ns = Bq and 3Bq. The mean 
energy of the resonantly scattered photons is reduced by drag when Fd mg /e\E\\ \ > 1 at z ~ \h. The critical temperature T& m „ above 
which this inequality is satisfied is plotted in Fig. [6] Although the drag formally is important when k^T^b > 0.2 keV, we find that 
other pair creation mechanisms will reduce the gap voltage in that temperature range and thereby reduce the effects of drag. 

The possibility of a surface gap in a magnetic field stronger than 4Bq deserve some comment. In this case most resonantly 
scattered X-rays will quickly convert to pairs, so that A' ~ 0. The gap structure examined here then appears to be unstable to 
collapse to small h, and it may disappear entirely. On the other hand, a self-consistent gap structure may be sustained outside the 
surface where B = 4-Bq, with properties very similar to those described here (Fig. 0J. 

3.3. Reduction of the Gap Voltage by 
Secondary Modes of Pair Creation 

The gap voltage can be reduced by additional modes of pair creation, which we now consider. Gamma rays created by the 
resonant scattering of thermal X-rays can collide with the X-rays to form free pairs, before they reach the threshold for single- 
photon pair creation (see Zhang & Qiao 1998 for a discussion of this process in the context of ordinary radio pulsars). Higher 
energy gamma rays will also be emitted even when B < 4Bq, with a cross section comparable to the Klein-Nishina value. The 
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FIG. 7. — Klein-Nishina cross section integrated over scattering angles that result in an outgoing photon above the threshold for direct pair creation, 
sinSj B2 ?i^ 7 > 2m e c 2 . The incident photon is assumed to propagate parallel to the magnetic field, which has only a modest effect on the cross section 

when the photon energy is larger than ~ 3TieB/m e c. Polarization effects are important only when the magnetic field is close to the threshold value 4Bq for direct 
pair creation following Landau de-excitation. In that case, an E-mode photon has a significantly higher threshold energy for pair creation than does an O-mode 
photon. The threshold condition for pair creation is essentially polarization independent when B <C 4Bq. 



emission rate is lower than for resonant scattering, but the conversion to a pair is almost immediate. This second process turns 
out to impose the stronger limitations on the thermal X-ray flux at low temperatures, and we discuss it first. 
The rate of pair creation integrated through the gap must not exceed a critical value, 



, hdpci _ A ( h 
dzn± < -—j— = — - — pa, (62) 



i7n±. < 

4 dr 4 v^ns, 

which corresponds to a net production of negative charges equal to the jump d24l i in the ingoing charge density at the outer 
boundary of the gap. We focus here on a gap situated at the surface of the neutron star, where outgoing and ingoing charges make 
significant contributions to the integral (1621 . 

3.3.1. Direct Pair Creation through Non-resonant Scattering 

We now work out the kinematic constraints on direct pair creation by non-resonant scattering. The cross section for scattering 
of a gamma ray by an electron is well approximated by the Klein-Nishina expression when uj > 3eB/m e c (Gonthier et al. 2000). 
The rate of Klein-Nishina scattering of peak photons is given by (Blumenthal & Gould 1970) 

2q2 

r™ = (1.65)^-J*ln( 7 e bb )il (63) 

(assuming that the thermal photons are emitted in a single polarization mode). Relative to the rate of resonant scattering of 
low-energy photons from the Rayleigh-Jeans tail, this is 

~ 2 x 10 -3 4t - 11 f[> ' Q (7 > 7®). (64) 

T- \i*J ln( 7 / 7 |) 17 1>} 

This expression holds in sub-QED magnetic fields for 7 > 7®(B/Bq) _i . 

In the rest frame of the electron, the X-ray photon has an initial energy Tiui' x ~ 7(1 - cos f^^Tiwx and is taken to scatter through 
an angle 6' kB 2 with respect to B. The threshold condition for direct pair creation following the scattering is 

Tilo'x 2 sin 6' kB 2 > 2m e c 2 , (65) 
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where 



W *' 2 l + (nuj' x /m e c 2 )(l-cosd' kB2 ) 
is the frequency of the photon after scattering. The quantity on the left-hand side of eq. ((65) is maximized when 

COS0' kB2 =— 1 2 , , ■ 

1 +m e cr /nw x 

The minimum energy of the incident photon for direct pair creation is 

huj x = (4+V26)m e c 2 . 

Approximating 1 -cos^b ~ 5<//^ns)~ 2 in the frame of the star gives the threshold energy 



7 = 



4(2 + v / 5)m e c 2 / r 



Tiuj x 



-V 



(66) 



(67) 



(68) 



(69) 



for direct pair creation by a relativistic electron moving away from the star. 

More energetic photons will scatter into a state that can convert directly to an electron positron pair over some range of 
scattering angles, 

Ax 2 [5 " Ax 2 [5 ' 

where x = 1 +m e c 2 /Tilu' x . The Klein-Nishina cross section integrated over this range of angles is plotted as a function of incident 
photon energy in Fig. [7] One sees that this integrated cross section peaks over the range of frequencies hu>' x = (10-20)m e c 2 . 
Denoting it by er n _ res , we make the approximation that the scattering is concentrated at a given frequency ftujQ ~ I6m e c 2 , so that 
o- n - res (ui x ) = a^_ le . & ui x 8 (uj' x - u>o) ■ Here trJJ_ res = 0.04 aj. (It will turn out that the rate of direct pair creation off a thermal photon 
field, maximized with respect to 7, is independent of u>q.) 
The rate of direct pair creation can then be written as 



2wdfl -^-{l T M) CTn-resC^X, M) 7)- 

hu>x 



(71) 



Here ujx and I UJx are the frequency and intensity of the thermal X-rays, which we assume to be emitted from the neutron star 
surface in a single polarization mode (the extraordinary mode), 

-1 



Tioji 



(2tt) 3 c 2 



exp 



1 



We evaluate this expression for values of 7 where u> = cJq/7(1 T A*) s i ts m tne Wien tail of the thermal X-ray distribution. 



x n-resv > U 



ft \47T 2 C 2 7 2 



exp 



ftuo/k B T hb 

7( 1 -Mmin) 



for an outgoing charge, and 



n— res 



(^,7) = 



ksTbb ( w o°n-re 



exp 



ftu>o/k B T hb 
2 7 



-exp 



ftujQ/k B T bh 

7 (1+Mmin) 



V 47T 2 C 2 7 2 / 

for an ingoing charge. 

The number of pairs that are created within the gap by each charge is obtained by evaluating the expression 



N± 



dz 
c 



(72) 



(73) 



(74) 



(75) 



where the Lorentz factor at height z is given by eq. d43l . The results are displayed in Fig. |8] where we have removed a factor 
h/R m - 

3.3.2. Collisions between X-rays and Gamma-rays 

The rate of pair creation by photon collisions can be written in terms of an integral over the distributions of the gamma rays 
and the secondary target X-rays, 



= I du> 1 I dojx / 2irdjjLx 



( dM~ 



\ du>^ 



/>«GJ ) ( ^-7- ) er X7 (w 7 ,Wx,/ix)(l-Mx)- 



(76) 



Here lo 1 = E 1 jft and ujx are the frequency of the gamma ray and the target X-rays. This expression applies to the case where the 
gamma rays are emitted parallel to B by outgoing relativistic charges. A similar expression applies to the case where the gamma 
rays are emitted by ingoing relativistic charges: one replaces /> with /<, and I- fix with 1 + fix- Here 
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FIG. 8. — Number of pairs that are created by the direct conversion of non-resonantly scattered photons, within the gap given by eq. i22\ . The background 
magnetic field is weaker than the critical value 4Bq for direct pair creation by reso nantly scattered photons. The contributions from outgoing and ingoing charges 
are evaluated seperately. The gap voltage is reduced by this process when eq. j87t is satisfied. Each set of three curves corresponds to magnetic fields 0.3Bq, Bq 
and 3Bq (thicker curves representing stronger fields). We choose parameters crj_ res = 0.04or, Hloq = 16m e c 2 , A = 1, Kims = 10 km, Rc = Rms and P = 6 s. 



The multiplicity dM 1 /dE 1 of resonantly scattered gamma rays is given respectively by eqs. (l58l and ( f60b in these two cases. 

We make the approximation that the cross section for collisions between X-rays and curvature gamma rays is concentrated at 
the threshold qjuix(1 - Hx) = 2m%, namely 

r /.„„2\2i 



m e c 



(78) 



The value of a Xl depends on the shape of the two photon spectra: following eq. (41) of Thompson (2008), we normalize to 
ox 7 = 0.02 oy (corresponding to a magnetic field B = 4-Bq). The integrals over uix, /ttx in eq. d76i l are then easily completed, 
giving 



= f>n G ] — ^—3 



cr x ^Q hb c f dM 



duj~ 



In 



1- 



exp 



2m e c 
huj y Q bh 



du>~, 



where we have taken /i m ; n = 0. The corresponding expression for the ingoing charges is 



: f<nci 



O"x 7 bbC f dM 



duj-, 



In 



1 - exp (-2m e c 2 //ju; 7 8bb) 
1 - exp {-m e c 2 /HuijObb) 



du>~. 



(79) 



(80) 



To obtain the net pair creation rate per unit area of the gap, the multiplicity dM. 1 /dE^ of gamma rays created by resonant 
scattering must be weighted by the distance Az = h — z (z) between the position z of resonant scattering and the outer (inner) 
boundary of the gap. The number of pairs that are spawned within the gap by each outgoing or ingoing charge is, then, 



1 



>.< />,<]»GJC 



doj~(h— Az) 



dri 



± >,< 



duj~, 



(81) 



Here Az is the distance from the injection boundary. The relation between frequency and position within the gap is obtained by 
assuming that Tvjj 1 = / reC oii7( Az)w e c 2 , where 7(Az) is given by eq. d43b . and the recoil factor /recoil by e q- (ED. The energy Tllox.q 
of the initial target X-ray is given by 



CJ > = 



fiux.o _ B/Bq 

^B?bb 7©bb 
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/recoilC-B/Bg) / ftuj. 



6bb \m c c- 

This expression applies to an outgoing charge; similarly one has 



1„ 



/recoil(fi/fiQ) / fiw. 



29 



bb 



(82) 



(83) 



for an ingoing charge. The integrals over lo 1 in eq. d79l ) and ( f80b can then be transformed to integrals over ui^ and w<. One finds, 



where 



_ 0.15a / eg b / _B_ 

^7 — i -r/m /recoil _ 



K; 



X 7 



'>,< 



(84) 



2 | "x 7 ^NS^C 

At 4 



x 7 



A 1/2 fL>zB 5 £ u \0.02a T ) V0.5keV 



* B 7b b \ 5 / P x 1/2 



and 



/[>= / rf£ >7 _^_in(l- e -">)ln 



6 s 



•Rc\ 



-1/2 



(85) 



' (Az/A) " 



duj^ , . ^ , : In 



' (Az/h) ' 



\- e -2"< 
\- e -u< 



l _ e -2u > /f, ccoil (B/B Q ) 

l_ e -42<l//r=»il(B/BQ) 
1,1 l_ e -25<,// rc cc,il(B/B Q ) 



(86) 



The minimum values of w> and u;< are found by substituting 7 = 7 max in eqs. (1821 . (1831 . The numbers of pairs that are spawned 
within the gap by outgoing and ingoing charges are plotted separately in Fig. [9] Note that we have extracted a factor of h/R^s in 
the definition of Kx 7 - 

3.3.3. Reduction in the Gap Voltage 

We now evaluate the effect of these secondary modes of pair creation on the gap structure. The voltage drops below the value 
given by eqs. (t40b and (l52l when the pair creation rate inside the gap is high enough that the inequality (l62l is violated. Self- 
consistent values of the gap thickness and voltage can be obtained by summing the rates of pair creation due to nonresonant 
scattering (j j3.3. It and collisions of gamma-rays with thermal X-rays (§ 13.3.21 ), and then arbitrarily adjusting $ gap by a factor £ gap . 
(The effect of resonant drag on the particle dynamics is neglected in calculating the pair creation rate.) A simple iteration gives 
a unique value of e gap at which the screening condition (l62l is marginally satisfied. The Lorentz factor at a fixed fractional depth 

z/h in the gap scales as 7 oc e gap , and the thickness as h oc e gap (eq. fl22l ). The coefficient Kx-y as defined in eq. (1531 scales as 

K Xl oc e gap . 

When the gap is partially screened, its voltage depends on the the outgoing and ingoing particle densities />.< x «gj (eq- 11771 ). 
The screening condition (|62t is readily expressed in terms of the quantities N± | and />,<, 



/>jv±l>+./W±l 

h/Rvs 



< 



(87) 



At the surface of the star, pair creation is mainly due to the ingoing charges when f < = 0(1). One expects that the gap voltage is 
only gradually reduced as the reverse inequality becomes satisfied: as the energy of the target X-rays is pushed further out onto 
the Wien tail, the number of targets drops exponentially. 

This expectation is borne out in Fig. [10] There is a critical surface blackbody temperature of about 0.15-0.2 keV above which 
the gap voltage is reduced. In the case where Z?ns = 3Bq, there is a reduction of a factor 3 at ^B?bb = -23 keV when J = 1.5,ogjc, 
and at k^T^ = 0.26 keV when J = \ .\pojc. A comparison of Figs. [8jand|9]shows that photon collisions are somewhat less efficient 
at making pairs than non-resonant scattering when Z?ns is stronger than 10 14 G, but are a much more important source of pairs 
when Z?ns is close to Bq and ^Ibb is larger than ~ 0.2 keV. The relative rates of pair creation by these two processes are plotted 
for a few values of Z?ns in Fig. QT| 

It should be noted that the total kinetic energy that is gained by the outgoing charges in crossing the gap, as given by eq. d52l 
with the correction of Fig. [TO] is larger than the power radiated in gamma rays during their transit through the gap. As a result, 
most of the pair creation that is induced by the gap occurs outside the gap, where the primary charges are decelerated by resonant 
drag (ED. 



18 



THOMPSON 



0) 
60 
t-i 
CO 

j= 

o 

0) 

a 

OJ 

CO 
0) 

l« 

o 



CO 



M 
O 



outgoing charge 


1 1 1 1 1 1 1 1 

^- - 


1 1 


ingoing charge 




— 












- 




**^<^**'^ 
* ^*>^ 










' ' - ' 




_ 


'/ ' 






/ / s 

' / ' 






/ // / 


B = (0.3, 1, 3)B q 




" ' / / 
<J / / / 
/ / / 
/ / / 
/ / / 


[light-heavy] 




~l / / 
/ / / 

' // / 













0.1 



0.2 



0.3 
k B T/keV 



0.4 



0.5 



FIG. 9. — Number of pairs created per charge flowing through a surface gap, by collisions between thermal X-rays and resonantly scattered X-rays. The gap 
solution is given by eg. 1221 . with B < 4Bq. The contributions from outgoing and ingoing charges are evaluated seperately. The gap voltage is reduced by this 
process where eq. (87) is satisfied. Each set of three curves corresponds to magnetic fields 0.3Bq, Bq and 3Bq (thicker curves representing stronger fields). We 
choose parameters A = 1, Rns = 10 km, Rc = i?NS and P = 6 s. 



4. CYCLOTRON DRAG ON THE SECONDARY PAIRS 
4.1. Increase in Pair Multiplicity due to Continued Resonant Scattering 

Pairs that are created close to the star will continue to scatter off ambient X-rays as they flow outward. The resonant drag 
force is stronger in more rapidly rotating neutron stars (e.g. young magnetars or young pulsars with high X-ray luminosities), and 
especially in magnetars with active, scattering magnetospheres. The flux of pairs can increase significantly through the conversion 
of the scattered X-rays, if the initial pairs are energetic enough. We neglect the possibility that the pairs are continuously 
reaccelerated, as they would be if the magnetic field were strongly turbulent (Lyutikov & Thompson 2005; Thompson 2008). The 
pair cascade then feeds off the initial kinetic energy of the outflowing charges, and additional generations of pairs can be created 
only within a limited distance from the star. 

Detailed expressions for the drag force experienced by relativistic magnetized electrons in a radiation field can be found 
in Dermer (1990), who was interested in the application to gamma ray bursts. Kardashev et al. (1984), Sturner (1995), and 
Lyubarskii & Petrova (1998) consider the interaction of outgoing particles with thermal photons near the polar cap of a neutron 
star, but do not consider the interaction further out in the magnetosphere, where the photon field is increasingly anisotropic. 

We suppose that outgoing electrons and positrons are created with 7^ > 7® (eq. l34l ) near the surface of the star. A charge 
will lose a significant fraction of its kinetic energy to resonant scattering over a distance ~ r if 



f pies 
/recoil 1 > 



1. 



(88) 



where / reco ii — B/Bq in magnetic fields weaker than Bq (eq. 11551 ). Evaluating expression for the scattering rate off a blackbody 
photon gas on the Wien tail of the distribution, one finds 



^=9.8x10^ f — 

ONS.14 \B, 



0.5 keV 



(89) 



■ #NS.14 V^Q, 

where w> =7>/7>- 

The resonantly scattered photons have an energy E 1 ~ 7 > (B/Z?Q)m ( ,c 2 . They continue to convert to electron-positron pairs off 
the magnetic field as long as 

E. ( 9 kB ~ E 1 — r — > 2m, c 2 . (90) 
Rc(r) 
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FIG. 10. — Reduction in the surface gap voltage 'I'gap due to pair creation inside the ga p by no n-res onant scattering of high energy X-rays, and collisions 
between gamma rays and thermal X-rays. The unscreened gap voltage & gll p0 is given by eqs. 440| and {52)- The result is given for various magnetic field strengths 
and for two values of J/pq]C. 



The pair creation zone sits close enough to the star that we assume that the field curvature radius Rc is dominated by the 
quadrupole/octopole component of the magnetic field (eq. HA9I ). Then the threshold condition d9"0l l can be expressed as 

3><#-f#V- (91) 



©bb \Bi^ t 

We can now determine the limiting radius R± for pair creation, and the characteristic energy of the outflowing particles at that 
radius. The value of lj^ deduced from eq. d89l varies from uj > ~ 10 near the surface of the star to £5> ~ 7 at the limiting radius 
for pair creation. The energy of the outgoing charges, as limited by drag, is 

7> (drag)= = B/Bq , (92) 

W>(1 -^minjibb 

and decreases from 7^ ~ 300(Z?Ns/3BQ)(A:B7bb/0.5 keV) -1 at the surface of the star. The magnetic field at r = R± works out to 

1/2 

s 0.5keV, 

Hence / reC oii — B^ m /Bq ^0.1, which means that the last generation of pairs that is created by the conversion of resonantly 
scattered photons is well separated in energy from the scattering charges. The characteristic energy of the latter has dropped to 

7,>(drag)=^ - 



B_(«i)=4.0x 10 I2 ^ 1/2 ( ) G (93, 



w > W >(1 -Mmin)©bb 
-1/6 R 2/3 

'-3 %S,14y 05keV 

and the secondary pairs have an energy 

7>[2] = ^=f?=)7>(drag) 



= 2 x 10 2 ^ 6 < 3 14 ( 7^Y /b (r = R±), (94) 



2m e c 2 \ 2B^ 

:10 ^ /3 < 3 -(oTk^)" 1/3 ir=R±y (95) 
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FIG. 1 1 . — Proportions of the pair creation inside a surfac e gap d ue to the direct conversion of non-resonantly scattered X-rays off the magnetic field ( j|3.3.U . 
and to collisions between gamma rays and thermal X-rays ( j|3.3.2t . The first process is more important at temperatures below 0.2 — 0.3 keV, and the second at 
higher temperatures. The current flowing through the gap is taken to be J = 1.5 Pgjc. The break in the curves at ^B^bb ~ 0.16 keV occurs where partial screening 
first becomes effective. 



The primary outflowing charges are injected with an energy 70 ~ 10 3 - 10 4 in the gap model described above. In this case, 
the outflowing charges multiply rapidly just outside the gap. Their energy spectrum can be derived as follows (again assuming 
negligible reacceleration). The mean number of photons scattered in a radial interval dr is dN 1 = Y^dr/c, and the loss of energy 
due to scattering is of 7 = (f KCO iij)T™ s dr/c. Focusing on the zone where B < Bq, and the gamma ray carries off a small fraction of 
the kinetic energy of the scattering charge, the energy spectrum of the gamma rays and created pairs can be written as 

dN± = dNj_ = ( t/[/ re coii7>( d ra g)] " 
'd^[2] 1 dE 7 \ d7[>(drag) 



7> [2] -7^=^-77^= , : ° • (96) 



Since 7 > (drag) scales approximately as ~ r , and B(r) oc r , this gives the asymptotic formula 

dN 4 

1^2}— (B«B Q ). (97) 
o?7[2] B/Bq 

Most of the pairs are created near the limiting radius R±, where the multiplicity of pairs per primary outgoing charge has increased 
to 

1 /3 

M ± „ 7o _ snn -2/3 -2/3 ( k B T bb \ / 70 \ 
M^Z " 2^ " 50063 BnS ' 14 lo^kevj ITo^J ■ (98) 

A beam instability will, in fact, be triggered in this outflowing plasma, since the distribution of particle kinetic energy extends 
to relatively low values. The effect on the particle energy spectrum will be strongest in the outer part of the pair creation zone, 
where / leco ii Cl. A calculation of the energy spectrum including the combined effects of resonant drag and the beam instability 
has never been carried out, but would be very useful. 

4.2. Limiting angle for Relativistic Pair Flows 

In this section, we make a brief detour to address the effects of resonant drag in two cases: i) the neutron star has a a relatively 
high X-ray luminosity (Lx > 10 34 ergs s" 1 ) and a faster spin than the known magnetars (P < 0.3 s); and ii) a significant fraction 
of the thermal X-rays are backscattered from the magnetosphere by cuiTent-caiTying charges. 

The strength and sign of the drag force acting on the outgoing charges depends on the angle 6^ between their motion along 
B, and the target photons. The thermal radiation emitted from the surface of the neutron star becomes increasingly collimated at 
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FIG. 12. — Close to the neutron star, the thermal X-rays emitted from its surface fill a wide cone. Beyond a distance R co i (eq. 11001 ), the angle between the 
propagation direction of a photon and the open magnetic field lines is dominated by the field-line curvature. 



greater distances; but at the same time the open dipolar 6 field lines become more strongly curved. The maximum angle value of 
QkB at the end of the open-field bundle is given by 

i_ M i max (%^y (99) 



2 V r 2 ' 4c J 

The first term is the large-radius expansion of eq. (|3T1 >. The collimation of the target photons about the local direction of the 
magnetic field is tightest at a radius (Fig. [T2b 



R - l= {mr s ) R ™ = UAR *i*{oT & ) R ™- (100) 

Expression d32b for the scattering rate r^ s of the outgoing charges is modified at r > R co i. One replaces 

< 1 -*-> 2 (101) 

The first factor represents the continuing dilution of the target radiation field as it flows away from the surface of the star, and the 
second represents the factor 1 - /i in the scattering rate. The fraction of the particle momentum lost to drag is given by 



/recoil T- S - = I ^ ) ( = (102) 

Note that the coefficient is independent of radius. If the spin of the star is fast enough, then the right-hand side of eq. ( 1102b can 
exceed unity - not only at r ~ R co \, but out to a very large distance from the star. The right-hand side is maximized at a}[> = 2, 
and is greater than unity if 

, / too \ 3/4 

,1/2 ( tx 

10 35 ergs s 



6 



We consider the effect of higher multipoles only when evaluating the threshold energy for single-photon pair creation at r < R co \. 
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This defines the critical spin period below which the motion of the secondary pairs is strongly reduced by resonant scattering at 

9 ^ $open- 

On a closed field line that is anchored at a polar angle 9(Rns) > (^ns/c) 1//2 , expression (11021 i generalizes to 

/reconr^ = ^&6\R NS ) (^j ule-*<>. (104) 

Charges moving through the closed magnetosphere will experience strong drag if the guiding magnetic field lines are anchored 
at an angle larger than 

/ too \ -3/8 

flmaxCRNs) = 1.36^S6 77^ \ de g" ( 105 ) 

■° V 10 ergs s 1 J 

4.2.1. Young Magnetars and Neutron Stars with 
Fast Spins and Large X-ray Luminosities 

Let us now consider the limiting Lorentz factor of the outgoing secondary pairs when cyclotron drag is strong at r > R co \. We 
focus here on the open magnetic field lines, and therefore on stars which are spinning faster than the known magnetars, with 
periods smaller than ( 11031 ). The force on a scattering charge vanishes when its Lorentz factor 7=1/ sin^g. The equilibrium 
value is 



1/2 



1/3 / \ -1/2 



r 



=12-4^j (r>* col ) (106) 

when the cyclotron scattering rate is very high. We can compare eq. ( 1106b with the Lorentz factor which the scattering rate of an 
outgoing charge is maximum, 



7^ _ B/B Q 



2 2(i-/Lt min )e b b 



1/3 / L oo \ -1/4 / r \-4 



=2.2xio^ N , 14 < 3 ,^j {^^—r) {^) ('>*-). 

(107) 

One has 7^, > 7 eq at r = R co i, so that resonant scattering decelerates an outgoing charge at this radius. The equilibrium Lorentz 
factor is attained a bit further from the star, at the radius 



R c ,= 23B ^^[0Al) {W^l (1 ° 8) 



where 

3/7 / f oo \ 1/28 

_ _ o , H-l/7 p-3/7 / _T \ ( Lx 



A charge is able to maintain 7 ~ 7 eq only over a modest range in radius. It continues to scatter off the Rayleigh- Jeans tail of the 
spectrum, and then decouples from the radiation field where the resonant frequency is uj^ < 1, and the asymptotic Lorentz factor 
is 

7oo = f -j-j 7eq(^eq) ~ (0.7-0.9)7 eq (^ eq ). (110) 

We reach an interesting conclusion: the outflowing charges are decelerated to a small Lorentz factor 7^ at a distance 20- 30/?ns 
from the star when the spin period is shorter than the value ( 11031 >. Unless the secondary charges are reaccelerated farther out 
in the magnetosphere, the radio emission will experience a high optical depth to induced scattering; and there would be a bias 
against detecting pulsed radio emission from neutron stars with high X-ray luminosities and short spin periods. 

4.2.2. Effects of Photons Scattered in the Closed Magnetosphere of a Magnetar 

The scattering of thermal X-ray photons by trans-relativistic charges in the closed magnetosphere creates a second component 
of the radiation field, which is weaker but much more isotropic (Thompson et al. 2002). Photons from near the thermal peak are 
scattered at a radius 

^(^) l/3 =^j n ^y 1/3 Rns . an) 

\m e cujx ) \ keV / 

This scattered radiation field can exert a strong drag force on outgoing charges. We derive an upper bound on its intensity, if the 
outflowing secondary charges are to avoid being decelerated to trans-relativistic speeds. 

The scattered X-rays subtend a wide range of angles, which means that they can resonate with the outflowing charges at a 
relatively small radius. Most of the drag force is supplied by photons from the black body peak, and so we approximate the 
spectral intensity of the scattered radiation as 

4 ^ 7 ^(¥)4^k)^-^ (112) 
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where Wbb ~ 1 keV in the case of an X-ray bright magnetar. 
The drag force in the rest frame of a charge moving parallel to B is 

F( { = J 'du'lit J ^a^i^dfi', (113) 

where the scattering cross section at the cyclotron resonance is given by 7 eq. (f30b . Transforming to the frame of the star gives 
fi'Iu' d\j! = j(}i—j3 e )I u! diJ,, and 

F,| = F ( = 7 T / d ^ ~ • (H4) 

Substituting expression (11 12b into this integral gives the drag time ?d rag = jm e c/\F« \, 

cfdrag(r) _( 1 2 \ 64u;l b Rl s (u; bb )m 2 e c 2 



f V Tres / 3cr7-B 2 ixr 

The scattered photons which resonate with the newly created pairs propagate in the direction 

eB 

l-/i= , (116) 

jm e cui bh 

and there is a maximum magnetic field for which a resonant interaction is possible (corresponding to /i ~ 0), 



B ma x = 8.6x 10 



2 (Z)(^) G. (117) 



V10 2 / VkeV. 
The drag force d 1 15b at B < Z? max is 

ctdngjr) _ Q 012r -i B lili5 R ^fi / 7 \ '/ 3 / ^bb\ 2/3 ^ B N _5/3 



res 



VlO 2 / VkeV/ Vfi ; 



(118) 



and depends weakly on 7. The resonating photons comprise a fraction ~ j(l - /j,) of the backscattered flux. The secondary pairs 
slow down unless r res Lx ^ 10 33 ergs s" 1 . 

5. LOW- VOLTAGE CIRCUIT WITH PAIRS SUPPLIED BY RESONANT SCATTERING 

We now turn to the case where the polar magnetic field of the neutron star is stronger than 4-Bq = 1.8 x 10 14 G. Close to the 
surface, a significant fraction of resonantly scattered photons are above threshold for conversion to an electron-positron pair. The 
distance Az' for the conversion of a gamma ray to a pair is therefore very small. The net effect on the surface gap structure is to 
take Rc — > in eq. d46i l. so that its thickness (|5TT i and voltage ( 152b are forced to zero. 

We therefore consider the alternative circuit model, with a lower voltage, that was outlined in §2.11 In this model, pair creation 
is smoothly distributed throughout the circuit, in such a way as to compensate the gradient in corotation charge density. The 
resonant scattering of X-rays from the exponential tail of the blackbody distribution can easily provide enough pairs: if the 
charges flowing through the circuit were energetic enough to scatter thermal X-rays from the black body peak, then the mean free 
path to scattering would be be very small, c/T Ks r <~ 10~ 2 - 10~ 3 from eq. fl35l >. 

The pair creation rate per unit volume is 

ri±~n > r™+n < r™ (119) 

in the zone where B > 4Bq. As before, and «<, are the number densities of ingoing and outgoing charges, and r'^ s , T r ^ s are 
the respective rates of resonant scattering as given by eqs. d32i l. d36i l. Substituting for ri± from eq. (fT3l . we find 

r-^+r-^^^i. (120) 

«gj «gj 2 dr 

We further have oc e" w> /7 2 > and r^ s oc e~ u}< /7<. The characteristic dimensionless frequencies £> and S< of the target 
photons are given by eqs. d33l and d38l . and are related by 

~ 1 — ^min ~ / 1 o 1 \ 

W > = UJ < . (121) 

In other words, an outgoing charge must attain a higher energy to scatter a photon of a given frequency than does an ingoing 
charge. This difference in resonant energy grows with radius. For a target photon of energy Tiux = ^B^bb, it is 



versus 

7^ = ^5xl0 2 fAWJ^r (123 ) 

v< 2e bb \Bq) lo.skevi ' 



7 Here we assume that the charge is moving relativistically, so that the photons are nearly collimated with B in its rest frame. 
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at r > Rns- 

Only a fraction of the outgoing and ingoing charges must scatter in order to satisfy eq. ( 1120b . which means that resonant drag 
has a small effect on the particle energies in this circuit model. The outgoing and ingoing charges pass through comparable 
potential drops and therefore have comparable energies, 7t> ~ 7<] . We deduce that T" s ^> T" s because ui^ > 2uj < and the 



(124) 



scattering occurs on the exponential tail of the Planckian distribution. Equation (1120b then simplifies to 



dig 
dr 



The right-hand side of this equation can be written as da dln(pa / B) / dr = AIqic / r, where the coefficient A (eq. 1201 ) includes 
the effects of relativistic frame dragging and magnetic field line curvature. The current carried by the ingoing charges in this 
circuit solution is given by eq. ( TTol l. The number of photons scattered per logarithm of radius works out to 



F^r dl a /dlnr 



I-Igj 



= A 



I 



(125) 



One of the charges in each created pair will typically be decelerated by the electric field and reverse direction before it can 
scatter another X-ray. Our toy model for the circuit assumes that this reversal is instantaneous, i.e. it occurs within a small 
distance of the point of creation. 

Substituting equation (|36l ) for F^ s in equation ( 1125l l, and taking into account that uj^ ^> 1, i.e. the resonance X-rays are in the 
far exponential tail of the Planck spectrum, we get the equation for 7<] , 



1 a f e bh r f B \ f B/Bq 

7 2 2% \B Q ) exp { 2 7< 6 bb 



= A 



I 



(126) 



For typical magnetar temperatures, k^T^ ~ 0.5 keV (8bb ~ 10 3 ), the solution of equation ( 1126b is well approximated by 



7< ~ 800B 15 



^B?bb 

0.5 keV 



(127) 



This corresponds to uj < ~ 15, and the scattering occurs in the exponential tail of the black body distribution, as expected. As the 
ingoing particles move to smaller radii, 7< increases in proportion to B oc r~ 3 . This acceleration is sustained by an electric field 
En . The equation of motion for ingoing particles reads 



2 d l< jrt 

-m e c = eh\\ 

dr 11 



Then we find 



eE» 



3 7< m e c 2 



The voltage distribution is found by integrating equation (11291 l, 



-e$(r) 



eEudr ~ 800m f c Bns.is 



0.5 keV 



Rns J 



(128) 
(129) 

(130) 



The net potential drop between the surface of the neutron star and the surface B = 4Bq is of the order of IBns.15 GeV. If the 
primary outgoing particles are electrons (pcj < 0), they are accelerated to the Lorentz factor, 



e$(r) 



(131) 



The outgoing particles undergo negligible scattering, and contribute little to the net rate of pair creation, since they see target 
photons that are a factor ~ 2 higher in frequency than do the ingoing charges. We therefore expect that this circuit solution 
should apply in the inner zone for either sign of pcj- 

The voltage ( 1130b can be compared with the gap solution ( 1521 , as reduced by secondary pair creation processes (Fig. |T0b . 
Normalizing the gap voltage to the threshold magnetic field Z?ns = 4Bq, we find 



$(r»fl NS ) 

$gapOBNS=4B Q ) 



4Bq 



£3 
0.3 



1/2 



P 

6s 



1/4 



0.05 (fe B 7bb = 0.1keV); 
0.1 0fc„7bb = 0.4 keV). 



(132) 



Seeds for the ingoing e ± discharge are supplied from the zone where B < 4Bq, and where voltage fluctuations are likely to be 
larger than they are in the inner circuit. In fact, the true solution in both zones may be time-dependent, with repeating avalanches 
of pair creation, as is seen in a one-dimensional model of the closed corona (Beloborodov & Thompson 2007). The main features 
of the steady circuit model should, however, survive in the time-dependent case: the ingoing particles dominate the rate of pair 
creation in the inner zone, that is, the pair avalanche is directed toward the star, and the time-averaged voltage is comparable to 
lGeV. 
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5.1. Stability of Voltage Solution I 

We explained in the preceding section why the surface gap described in $3]cannot be maintained when Z?ns > 4Bq. We now 
address the stability question from the opposing angle: when Z?ns < 4Bq, does a gap spontaneously form in a circuit that is set up 
in the low-voltage solution, with smoothly distributed pair creation? We approach this question here by considering the stability 
of electrostatic modes in the circuit, but note that other degrees of freedom could be involved. In particular, a high-frequency 
current fluctuation would be generated by a turbulent cascade if the magnetic twist received a low frequency perturbation of a 
moderate amplitude (Lyutikov & Thompson 2005; Thompson 2008). 

The instability in question bears some resemblance to the formation of a double layer within a current-carrying plasma. In 
an electron-ion plasma, such a structure is believed to be seeded by a two-stream instability, such as an ion-acoustic instability, 
which is driven by the drift motion of the electrons that supplies the current. Localized wavepackets that are formed via this 
instability undergo non-linear growth, leading to the formation of plasma holes which are bounded by layers of positive and 
negative charge. Multiple double layers can form within a single circuit, as is observed in the Earth's magnetosphere. A diode 
structure with a deeper potential drop may form through the coalescence a large number of weak double layers. 

The low-voltage circuit consists of counterstreaming clouds of electrons and positrons, and so we first review what is known 
about the non-linear evolution of the two stream instability in a strongly magnetized pair plasma. Near the surface of a neutron 
star, the electron cyclotron frequency is many orders of magnitude larger than the plasma frequency, and so filamentation insta- 
bilities are suppressed. Guidance about the behavior of the two-stream instability can therefore be provided by one-dimensional 
simulations of the purely electrostatic instability. The perturbation will, in practice, have an electromagnetic component because 
plasma oscillations on different field lines have different phases. 

We imagine preparing the ingoing and outgoing charges with a distribution function concentrated at momenta pu = ±jVm e 
parallel to B. For simplicity, we focus on the case where pcj — > (e.g. the open flux tube is anchored near the rotational equator). 
Then the plasma frequencies of the electrons and positrons are equal, and are given by 

2 4717V 47!7V 2irJe 
uj p = = = , (133) 

m e m e m e c 

where the densities refer to the frame of the star. The dispersion relation is obtained from the longitudinal part of the dispersion 
tensor, 

1- = - , - . - =0. (134) 

The fastest growing mode has a growth rate and wavenumber 

r. = ^j; *. = 3*£, (135) 

and a vanishing real frequency. Purely growing modes exist within the range of wavenumbers < k < 2 3 / 2 r*/V = 1.6£». In 
practice, there is an imbalance in the densities of positive and negative charges, so that the growing modes have real frequencies 
of the order of T* and finite phase speeds V P h ~ c. 

Simulations of colliding relativistic plasma clouds in one dimension 8 (e.g. Dieckmann 2005, Dieckmann et al. 2006) show the 
growth of electrostatic modes with a modest range of k centered on which then develop into a broad spectrum of electrostatic 
turbulence on a timescale ~ 50 T~ l . A necessary condition for the existence of a beam-driven mode is, therefore, 

50^</? NS , (136) 



or equivalently 



I 



,< (l0-^) 2/3 =7xI0XS„< 6 (£)"" 3 (£)' P - 

Comparing this expression with eq. ( I127l i for the equilibrium Lorentz factor 7<, = 80Z?NS,i4(^B7bb/0.5 keV)" 1 of the pair-creating 
charges (eq. 111271 ). one sees that this inequality is easily satisfied. Note that the growth of the beam-driven mode is not suppressed 
by the radial inhomogeneity of the pair plasma, because growing modes exist over a relatively broad range of wavenumbers. 

A similar inequality can be obtained by supposing that the distribution function is flat at small momenta, so that the average 
of 7 -3 over the distribution function is dominated by a small fraction ~ I/7 of particles with momenta pu ~ ±m e c, (7~ 3 ) 1//2 ~ 

1 /7 1 / 2 . A transrelativistic particle propaga ting t hrough the background electric field gains energ y at the rate ^dj/dt ~ 7c/J?ns> 
which when equated with the growth rate ( 11351 ) gives the same dependence on 7 as in eq. ( 11361 ). 

The rate of pair creation changes in response to the beam-driven mode, and screening of En by the pairs can suppress growth. 
We imagine that gamma rays are emitted and converted to free pairs, and that one of the particles so created immediately reverses 
direction in the background electric field. There is a critical value of the total free path for this process, denoted by Az, below 
which the screening effect is important, 

^3/2 / p \ 1/2 

Az < k~J ~ f/ 2 — = 0.5 -475— r-l cm - ( 138 > 

Up Tj l/Z \6 s / 

8 These simulations are of two colliding electron-ion clouds, in which the ions supply a positive neutralizing charge distribution and do not strongly influence 
the growth of the instability. 
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When this inequality is violated, the gamma rays that are emitted within a localized potential perturbation convert to pairs at 
distant points in the circuit, where they cannot influence the growth of the potential. The minimal rate of pair creation needed to 
screen the growing, beam-driven potential fluctuation 6$ before a broad spectrum of electrostatic turbulence develops is 

eSri ± > ^Jp- (139) 

Here 

k 2 S$ 

Sp=~f L , (140) 

4-TT 

is the corresponding charge density fluctuation. The change in the rate of gamma-ray emission per unit volume resulting from a 
potential fluctuation <5$ is given by 

dn\ e<5<i> 

<W 7 = 1 r^ J' < 141 > 

07 m e c A 

If inequality ( 11381 ) is satsified, then 8n± ~ <5n 7 , and we have a second condition for effective screening 



fc»<(10 2 — ■ ( 142 ) 
\ m e c 3 07 J 

The background rate of pair creation is given by eq. (1120b . namely eh± ~ (c/2)dpci/dr = (Ac/2R^s)pcs- Substituting for into 
eq. (11421) . and neglecting the difference between J and pac, one finds 

/ 1 u pRns \ 2/7 (dlnn ± yV 7 2/7 1/7 1/7 / P V 1/7 ( d\nh ± \ ^ 

The screening of the growing potential perturbation requires that both conditions ( 11381 ) and ( 1143b be satisfied. The second 
condition is generally satisfied by our circuit model. The first condition is easily satisfied if Z?ns > 4Bq, so that pair creation 
follows resonant scattering; but it is not satisfied if B < 4Bq. Then the two-stream instability in the pair plasma will be effectively 
triggered near the neutron star surface, and the non-linear growth of the potential fluctuation into a strong diode may result. 

6. DISCUSSION 

We now confront the circuit physics described in this paper with a few key observational tests, and then summarize our 
conclusions. The behavior of the radio magnetars appears to be unusual in several respects in comparison with long-period radio 
pulsars. 

6.1. Application to the Radio Magnetars XTE J1810-197 and IE 1547.0-5408 

6.1.1. Radio Luminosity 

The spindown of XTE J18 10— 197 is transiently accelerated (Camilo et al. 2007a). The polar magnetic field that is deduced 
from the standard magnetic dipole formula has decreased from Z? po i e = 2Z?mdr ~5x 10 14 G to 3 x 10 14 G, and a further decrease 
is possible. The polar field of IE 1547.0-5408 is inferred to be Z? po i e ~2x 10 14 G. As a result, the actual polar fields of these 
neutron stars could be weaker than 4Bq, in which case the gap model described in $3]is relevant. We first consider whether the 
power dissipated in the gap can supply the observed radio output. 

When the detected pulsed radio emission of XTE J 1 8 1 0- 1 97 was at its peak, the Goldreich- Julian charge flow on the open field 
lines was Ics/e ~ BMDR^NS^ 2 / ec = 2 x 10 31 s" 1 . The mean energy per outflowing particle crossing the gap is e$ gdp = 7max'WfC 2 ~ 
5 x 10~ 3 erg (eq. l52l ). including a factor ~ 0.5 reduction in the gap voltage due to the screening effects of secondary modes of 
pair creation at a temperature ^B?bb ~ 0.25 keV (Fig. ITOb . The net dissipated power is therefore Ica&gap = 1 x 10 29 ergs s" 1 . The 
peak isotropic output at 40 GHz is, by contrast 4nD 2 vS u = 1 x 10 30 ergs s" 1 at v = 100 GHz at a distance D = 3 kpc (Camilo et al. 
2007a). The apparent radio output could be enhanced by an order of magnitude by the effects of beaming, but even in that case 
the two results are barely consistent. 

We must therefore consider two alternative possibilities: that the same gap-like structure can form on a broader bundle of 
magnetic field lines, extending into the closed magnetosphere; or that the power dissipated on the open field lines is much higher 
than is implied by the gap model of $3] It should be kept in mind that the spindown power of XTE J18 10-197 was at least 
2.4 x 10 33 /ns 45 ergs s" 1 at its peak, and so it can easily accommodate the radio and infrared output (Lir ~ 1 x 10 31 ergs s" 1 : Rea 
et al. 2004; Wang et al. 2007). 

6.1.2. Peak Emission Frequency and Pulse Width 

The pulse duty cycle of the two radio magnetars is much larger than would be expected if the radio emission originated close 
to the star on the open magnetic field lines. We deduce either that the emission zone is close to the light cylinder, or that there is 
coherent emission from the closed magnetic field lines. In the case of IE 1547.0-5408, this effect could be partly explained by a 
nearly alignment of the magnetic axis with the rotation axis (Camilo et al. 2008). Nonetheless, the fit to the linear polarization 
curve found by these authors does not require that the tilt angle between fi and fi be nearly as small as 6 open (R^s), if the radio 
beam is intrinsically broad. (Indeed, the small apriori probability of detecting a radio beam of angular width ~ 26' op e n (^Ns) from 
this object can be used to argue that the beam must be broad.) 
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It is possible to deduce a value of -fA4± under the assumptions that i) the radio waves are emitted tangent to the magnetic field 
lines within a maximum angle which is a multiple Ne of 9 open {r); and ii) the emission frequency at radius r is equal to a multiple 
N y of Vp±, We consider an outflowing gas of pairs, within which the plasma frequency is 



4 ± = fM^I (144) 

Ztt \ irm e J 

in the bulk frame. The proper density of electrons and positrons is 

/ 1 M± 
2n± = -{n e + + n e -)= n G j (145) 

7 7 

where Ai± ^S> 1 is the number of light charges created per primary Goldreich-Julian charge. We assume for simplicity that the 
plasma is in bulk motion with a Lorentz factor 7, so that the frequency is increased in the stellar frame by a factor 2j, 

, o„ ± , 2 N 1/2 



' '2nLe^ 1/2 



— ■ (146) 



The plasma frequency is maximized close to the surface of the star, 

1/2 

^ e i(flNs)=' 



2ileB NS 



n 2 m e c 



= 2.6 xlO 17 - — ±\ — — Hz. (147) 

V 10 4 J \AB Q J \6sJ 

Outside the maximum radius of pair creation R±, one obtains the usual scaling, 

/ \ "3/2 

^±W=f^J v?±(R±) (r>R±). (148) 

Note that depends on the product jAi±, which is approximately invariant if the charges undergo resonant scattering and the 
scattered photons are able to re-convert to pairs. In the gap model described in $3] one has jM± ~ e$ ga p/m e c 2 . 

The tangent to a dipole field lines makes an angle ~6 with respect to the magnetic axis when 9 -C 1, and so the pulse duty cycle 
St/P is given by 

St = 2Ng (§6>op en ) =N 36 opS n 
P 2tt 2tt 

The multiplicity corresponding to a duty cycle St/P at a frequency v = NuVp]. is 

7A1±= W' ch '1"stJ [r.) [as-J ■ (150) 

Substituting numbers appropriate to XTE J 1 8 1 0- 1 97 (St jP = 0.2 at v = 2 GHz and P = 5 .44 s) gives 

7 X± -5 x 10 8 AT 2 A^ 6 (XTEJ1810-197). (151) 

Let us first compare this result with the turbulent flux tube model that is described in Thompson (2008). In that model, the pair 
multiplicity is jA4± ~ (elcs/m e c 3 ) = 1 x 10 8 . The plasma frequency at the base of the flux tube then sits in the optical-IR band, 
and it is possible for the optical-IR emission of XTE J18 10-197 to be an upward extension of the spectrum at ~ 100 GHz. If this 
were the case, a narrower pulse profile would be expected in the IR band as compared with the radio band. 

The alternative explanation for the broad pulses and high intensity of the radio emission of XTE J18 10-197 is that the outer 
closed magnetic field lines support a unidirectional discharge. Within a cr itical angle 9(R NS ) ^ 3°(L X /10 34 ergs s" 1 )" 1 / 2 ~ 
9# pen(^Ns) from the magnetic axis, outflowing relativistic charges are not stopped by resonant scattering of the thermal X-rays. 

The gap model of Ogives r )M.± ~ 1 x 10 4 , which is consistent with eq. ( 1151b if N^ 2 Ng ~ 6. The observed radio power is also 
reproduced, because the cross-sectional area of the emitting field lines is 30-100 times larger. 

Although we cannot in this paper decide between these two possibilities, the spectrum of the radio magnetars provides an 
important clue. The very high peak frequency of the emission does not, by itself, demand a high plasma density on the open 
magnetic field lines. But the comparison with ordinary radio pulsars (for which the peak frequency is ~ 0.1 — 1 GHz) does 
suggest that the plasma density is much higher in the radio magnetars than in pulsars, by a factor > 10 4 — 10 6 if one uses the 
scaling n± oc ^ 2 eak . 

6.2. Is 4Bq = 1.8 x 10 14 G a Limiting Polar Magnetic Field for Radio Pulsars? 

The distribution of neutron star magnetic fields, as inferred from pulsar spindown, offers important tests of the physics ex- 
plored in this paper. Considering only sources that are discovered in the radio band, 9 the distribution of dipole fields has an 

9 As noted above, the polar fields of the radio magnetars cannot be deduced with precision from their spindown rates due to the presence of large-amplitude 
torque variations, but are not obviously inconsistent with this bound. 
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upper envelope at Bmdr — 1 x 10 G. The corresponding polar magnetic field is remarkably close to 4Bq = 1.8 x 10 G. A 
previous suggestion that radio emission could be quenched in super-QED magnetic fields (Baring & Harding 1998) depends on 
assumptions about the polarization dependence of photon splitting that are not born out by microphysical calculations (Adler 
1971; Usov 2002). Our focus here is different, on the voltage of the open-field circuit and the properties of the gap that may form 
there. 

When the surface magnetic field is stronger than 4-Bq, the open-field voltage depends on whether a gap is able to form displaced 
from the surface of the star. The most plausible mechanism of gap formation involves the generation of multiple weak double 
layers in the circuit, seeded by a two-stream instability, which then accumulate into a diode structure with a larger voltage. The 
two-stream instability is stronger in the low-voltage circuit solution described in ^5]than it is the higher-voltage solutions with a 
gap, but the absence of a sharp boundary at the surface B = 4Bq may impede the second stage of the process. 

Whether or not a gap does form in this situation, pairs are easily created by resonant scattering throughout the inner circuit. 
Given that the net difference in circuit voltage between the two solutions is only a factor of 10, the detectability of radio pulsations 
from quiescent neutron stars with Z?ns > 4Z?q and P ~ 1 — 10 s is probably a matter of flux sensitivity, rather than the presence or 
absence of emission. 

6.3. Summary 

1 . - We have described a self-consistent circuit solution with a gap at the surface of the neutron star, which requires a current 
density \J\ > pcjc and a surface magnetic field Z?ns < 4Bq. The voltage across the gap is e|$ gap | ~ 10 4 m e c 2 when k^Tu, < 
0.15-0.2 keV. We have reconsidered the effects of magnetic field line curvature on the voltage, and find that the octopole 
component makes the largest contribution to linear order, with the quadrupole contributing to second order, when these two 
components are treated as a perturbation to the dipole. It is not known whether the two detected radio magnetars have polar 
magnetic fields stronger or weaker than the critical value 4Bq, after allowing for changes in magnetospheric structure resulting 
from the injection of currents during periods of X-ray activity. 

2. - The dominant pair creation mechanism within the gap is the conversion of resonantly scattered photons off the magnetic field 
when the blackbody temperature is relatively low, ^B^bb 0.2 keV. At higher temperatures, the voltage is reduced by secondary 
modes of pair creation. Up to a temperature of ~ 0.3 keV, the most important such mechanism involves the non-resonant (Klein- 
Nishina) scattering of X-rays, followed by direct conversion off the magnetic field. At yet higher temperatures, or in fields weaker 
than ~ 10 14 G, collisions between gamma rays and thermal X-rays occur at a higher rate. The net effect is to reduce the gap 
voltage by a factor ~ 0.1 when ^B^bb ~ 0.3 keV. 

3. - Screening of the electric field outside the gap can be effected by the polarization of the outflowing pairs. The pairs injected 
near the outer boundary of the gap continue to downscatter off the ambient thermal X-rays, thereby creating secondary pairs. We 
have deduced the spectrum of energies that result. Inside a distance r ~ 10/?ns, the outflowing pair cloud maintains a charge 
density equal to pa, as charges of the opposite sign reverse direction and return to the star. Mildly relativistic electrons and 
positrons feel a strong drag force near r ~ 10/?ns, but a steady, charge-balanced flow remains possible in which the 'negative' 
charges flow outward sub-relativistically. 

4. - Where B > 4Bq, the open-field circuit sustains a much lower voltage, because resonantly scattered photons convert directly to 
free pairs. In the temperature range k^T^ ~ 0.1-0.5 keV, the voltage is reduced by a factor ~ 0.05-0.1 (Bns/4Bq) as compared 
with the gap model evaluated at the threshold field Z?ns = 4Bq. We have also derived the voltage distribution through the inner 
circuit in the case where pair creation is smoothly distributed, at just the rate that is needed to compensate the gradient in pcj- 
There may be a connection between the quenching of a surface gap at Z?ns > 4Bq and the absence of radio pulsars with dipole 
fields stronger than 10 14 G. 

5. - The radio output of XTE J18 10-197 and IE 1547.0-5408 is somewhat brighter than is allowed by the gap model of §f3] if 
the emission is restricted to the open dipolar field lines. The radio pulses are also much broader than is typical of long-period 
radio pulsars (Camilo et al. 2006, 2007b). The same gap model predicts much narrower pulses than are observed if the radio 
emission frequency is close to the local plasma frequency, and the emission comes from the open magnetic field lines. When 
combined with the very hard radio spectra, this suggests that the rate of pair creation is much higher in the radio magnetars than 
in more typical radio pulsars. The wide pulses, hard spectra, and rapid variability of the radio magnetars are more consistent 
with the electrodynamic model described in Thompson (2008), in which a broad spectrum of Alfvenic turbulence is excited by 
current-driven instabilities in the outer magnetosphere. 

6. - Relativistic pairs flowing outward from an X-ray bright magnetar will suffer strong resonant drag beyond a radius R co \ ~ 
50(P/6 s) 1//3 i?Ns if the guiding field line is displaced too far from the magnetic axis. At a spin period of 6 seconds, this critical 
angle is ~ lO0 O pen; but if the spin is faster than ~ 0.4(Lx/10 34 ergs s" 1 ) 3 / 4 s, then most of the secondary pairs on the open field 
lines will suffer strong drag outside ~ 10 stellar radii. Radio emission from these charges will be suppressed by induced scattering 
unless they are continuously reaccelerated. There is therefore a bias against detecting pulsed radio emission from neutron stars 
with moderately fast spins and high X-ray luminosities (e.g. young magnetars). 

This work was supported by the NSERC of Canada. It began as an investigation of the radio death line of magnetars, and I am 
indebted to Andrei Beloborodov for his contribution to those early efforts. In particular, he provided a cogent argument that pair 
cascades driven by curvature radiation cannot be self-consistently sustained in magnetar magnetospheres. I also thank him for 
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APPENDIX 

SOURCES OF GRADIENT IN COROTATION CHARGE DENSITY 

The voltage on the open magnetic field lines is dominated by the gradient in the corotation charge density in the case where the 
first order difference between J and pcic is compensated by a returning flow of charges of the sign opposite to the local corotation 
charge. Here we compare the effect of general relativistic frame dragging (Muslimov & Tsygan 1992) with field line curvature 
as dominated by a quadrupolar and octopolar magnetic field. 

A purely dipolar magnetic field in a Schwarzschild spacetime has the shape 

B(r)=-^-B dipole (r), (Al) 

where f(r) = 1 + £^ j|g(2GMNs / 'c 2 r) k (Muslimov & Tsygan 1992, eq. 32). The prefactor here describes the deviation of the 
poloidal magnetic field lines from the shape they would have in the absence of gravity. It does not affect the voltage directly, 
because both the charged particle flux and the magnetic flux satisfy the same conservation equation. The relativistic charge 
density a{r)pci(r) has the form (0, taking into account the effects of frame dragging. When evaluating the lapse function, the 
stellar rotation can be neglected: on has 

/ 2GM Ns y /2 

a{r) =1 = (A2) 



outside the star. It is sufficient to take a — *■ 1 and / — * 1 in all expressions, because the magnetic flux density and the relativistic 
current density J = apy and the magnetic flux density satisfy the same conservation law 10 V ■ B = V • J = 0. 

An additional mechanism for voltage generation becomes significant for a more general magnetic configuration in which the 
open field lines are strongly curved. The ratio pa/B increases with distance from the star on those field lines that bend toward 
the spin axis. In the absence of pair creation this implies an increasing mismatch \p\ > |/Ogj,o|, and a voltage is generated that 
accelerates the Goldreich- Julian flow (e.g. Barnard & Arons 1982; Asseo & Khechinashvili 2002). This effect dominates over the 

10 Here V is the gradient operator in Minkowski space. 
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frame-dragging effect beyond a certain radius, and can be especially important in magnetars due to the screening of the parallel 
electric field by direct pair creation close to the star. 

We focus here on the case where the magnetic field is predominantly dipolar beyond zone I. The presence of a higher multipole 
(quadrupole or octopole) causes a bending of the open dipolar field lines. The angular deflection 86{r) can be deduced by 
working inward from the light cylinder. The radial component of the magnetic field is dominated by the dipole component, 
Bq = B^s(r / Rns)~ 3 ■ Well inside the light cylinder, but still far from the magnetar surface, the non-radial field component on the 
open-flux bundle is dominated by the higher multipole, e.g., 

/ r \- e ~ 2 ( r \" £ " 2 

B e (r) ~ B e (R NS ) I — I = F> Bxto I — I 



; £e Bns 



We suppose that this additional multipole is oriented randomly with respect to the magnetic dipole axis. Then 



(A3) 



S0(r) : 



-£+\ 



(A4) 



Bg(r) dr _ e t t r \ 

B D (r) r~l-\ \R m J 

The parameter ee contains the angular dependence of the £-pole. 

The lowest-order multipole component which makes a significant contribution to voltage and curvature of the open field lines 
is the octopole. 11 To see this, note that the transverse deflection of the open magnetic flux bundle away from a pure dipole is 
L± = r86(r) cx r 2 ~ l , which is constant for a quadrupole. More generally the curvature radius can be found from 



which in combination with 



gives 



Rc = 



B~ [\-h e (rf] 1 ' 2 f+b e (f)6; 



\(B-V)B\ 



be(r) : 



B D (r) 



d(rbe) . 



to second order xnbg. A single additional multipole imparts a curvature radius 



Rr = 



Rns 



(£-2)(e-l)S8(r) (l-2)et 
to lowest order in be- This expression gives infinite Rc for 1 = 2, and 



- 

Rns 



Rc 
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(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



when an octopole is added to the dipole. The quadrupole provides a second-order correction to the octopole curvature, which has 
the same scaling in radius, and is of comparable magnitude if £3 ~ e\. 
When the magnetic field is mainly dipolar, the opening angle is given approximately by the dipole formula, 



/open(x)" 



fir 



nl/2 



= 0.014« '(J-) 



-1/2 



1/2 



rad. 



(A10) 



The factor / open (x) = 0(1) corrects for an inclination angle x 7^0 between the magnetic moment and the rotation axis; we generally 
set / p en = 1 in expressions in this paper. The tilt of the field lines with respect to the magnetic dipole axis can easily be much 
larger than 8 open . For example, a quadrupole/octopole component of magnitude 



^3 : 



:£ 3 + e 2 ^>0.012/ 1 / e 2 n < 2 6 



P 

6s 



-1/2 



\Rns) 



1/2 



(All) 



is sufficient to impart a tilt SO > # open at a radius r. In this situation, both 66 and pes are approximately constant across the open 
flux bundle. 

The electrostatic potential, as given by eq. (1A181 >. depends on the deviation of the charge density from the local corotation 
density. 12 In this case, there is a direct relation between the curvature of the open flux bundle and the gradient in pa parallel to 
the magnetic field. One has 

d( Pa /B)_ n [(j g. V)j g ] = _^_ sinxcos0C; (A12) 



di 



2lTC 



2ttcR c 



11 To linear order in the multipole amplitude. This conclusion presupposes the dominance of the dipole component of the magnetic field in the acceleration 
zone. If the dipole component were absent, then a quadrupole by itself would still not supply significant field line curvature. The opening angle in that case would 



be O , 



' flr/c, and the curvature radius would be comparable to the light cylinder radius, Rc ~ >"/6> op cn ~ c/fi. A mixture of multipoles is required to provide 



an effective particle accelerator. 

12 For simplicity we set the lapse function a 



1 in what follows. 
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or, equivalently, 

d(p G! /B) Pgj . 

B = — sinxcos^c- (A13) 

dr R c 

We have made use of the fact that (B ■ V)B points in the ^-direction in the plane of curvature of the open magnetic field lines. 13 
Here 4>c is the angle by which the plane of curvature is rotated with respect to the ft plane. The zero of <fic is defined so that 
the field lines bend toward the rotation axis when cos <\>c > 0. There is no acceleration when cos <f>c < 0. 

Limiting Voltage on a Slender Magnetic Flux Tube. 

For completeness, we examine the voltage in an extended gap which has a length Ar much larger than the radius 9 open r of the 
open magnetic field bundle. When the star rotates steadily and the magnetic field is fixed in the corotating frame, an electrostatic 
potential $(r) can be introduced, which satisfies the equation, 

V • (a 1 V$) = -4tt(p- pa). (A14) 
Then the electric field component along B is given by 

£j| =oT 1 B- V$. (A15) 

Once again, we make the approximation a — ► 1. We approximate the open flux bundle as being locally cylindrical, with a 
cross-sectional area 

S^~ir9 2 opea r 2 . (A16) 
Choosing a longitudinal coordinate r and transverse (cylindrical) radius 03 = Or, Gauss' law can be expressed as 



IJL 

03 903 



d 2 <5> 
'd? 



+ —=-4tt(p- Pgs ). (A17) 



It is usual to apply the boundary condition $ = at 03 = 6* open r, which must however be reconsidered when the closed magnetic 
field lines are twisted by the release of internal stresses. In this situation, a modest voltage ~ 10 2 - 10 3 V develops on the closed 
field lines, as is required to sustains the current (Beloborodov & Thompson 2007). The power dissipated in the closed-field 
circuit is drawn from the energy of the toroidal magnetic field, and the gradual relaxation of this twist is associated with a time 
derivative of the longitudinal component of the vector potential. The twist on the open magnetic field lines is, by contrast, 
continually replenished by the rotation of the star, and decreases only adiabatically as the star spins down. We therefore retain the 
boundary condition $ = at 03 = 8 open r, but it should be understood that this boundary condition could be modified if the outer 
magnetosphere were subject to fast, current-driven instabilities. 

A simple solution to eq. JA17b is obtained by neglecting the longitudinal component of the laplacian and taking both p and pcj 
to be constant across the open flux bundle. This approximation is reasonably accurate outside a distance r-R^s ~ (6* pen/2)^NS 
from the neutron star surface. The solution is 




*(r,ff) = — [ 1 j ( 1-— i • (A1S) 

Here 

/ G j(r) = S ± (r)a(r)p C] (r)c (A19) 

is the net current along the open field-line bundle. The twist on the open magnetic field lines is supported by the outward flow of 
the corotation charge. Neglecting the effect of frame dragging, the charge flow rate is given by 

^ = ^^ns = 4x10 3i B ns , 15 4s,(^)" 2 a" 1 - (A20) 

An important feature of eq. (1A181 > is the reversal in sign compared with the solution that would be obtained by neglecting the 
transverse component of the laplacian in favor of the longitudinal component, and applying the boundary condition $(/?ns , 0) = 0. 
This means that if the charge density p decreases with radius more rapidly than pes, then an electric field is generated that is of 
the right sign to drive the charges away from the star (Scharlemann et al. 1978; Muslimov & Tsygan 1992). 

We suppose that screening is absent beyond some radius Enough charges are created inside this radius so that p = pes', but 
that the charge density is determined by the continuity of the outward charge flow farther out in the magnetosphere. The voltage 
solution beyond this radius is, then 



1° 

$(r,0) = K*-^ 



^ -1 



open 



Here we have used 



K , = K (R,)=^1. (A22) 



' To leading order in the quadrupole/octopole amplitudes; eq. IA7I . 
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Expanding equation (IA2 11 > near gives 



7° / Ar' 

+ Ar) -$(£*) ~ -3k* — — )■ (A23) 



c \R* 

The limiting Lorentz factor of an electron (or positron, if pes > 0) flowing beyond r = is 



7max(fl*) = ^ = — ^ 



GI NS n 2 eB(R*) 7 

= Z = 5.3 X 10'/ NS .45 



10 15 G 



P 

6s 



-2 



(A24) 



Now let us turn to the case where the gradient in pa is dominated the quadrupole and octopole magnetic field. The voltage 
that develops in the absence of screening can be obtained by integrating eq. (IA12b along B, 



-t+\ 



r 



r>/?„, £>3. (A25) 



In this expression, p = j^ns^ns * s tne ma g net ic moment of the star, and 

ne = {£-2)eg sinx cos</>c- (A26) 

The expansion near r = R* is given by 

1° / R \ ( Ar\ 
+ Ar) -$(*.) ~ K£ -21 f-JL j j, (A27) 

The limiting Lorentz factor of accelerated electrons (or positrons, if pcj > 0) that can be accelerated outside a radius R* is 



7max(P*) = . \ 3 I 1 ■ (A28) 

£- 1 m^c- 5 \R NS J 

We can now compare the voltages that are induced by the frame dragging and field-line curvature effects. Outside a radius 
> /?ns, one has 

$(field curvature) Kg ' " " " ' " ' 



P* \ g ^ ^£ / P* \ (A29) 

^?NS / ^-1^NS,45 \RnS J 



following eq. JA9I ). One then has 



$(frame dragging) (^ - 1)k(%s) 
In the case where quadrupole/octopole components are present, one can define an effective octopole amplitude 

~ / de 2 \ 

k 3 = I e^ + e 2 — j sinxcos0 c , (A30) 
$(field curvature) /?ns.6 / ^* \ 

= 3.4^3 . (A31) 

$(frame dragging) 7 NS , 4 5 \Rns J 

In the particular case where B = 4Bq, the quadrupole/octopolar components of the magnetic field will dominate the effects of 
frame dragging if k 3 > 0.17 B^ 3 15 I nSt 4 5 R^ s 6 . 

A deviation of p from pes also arises due to the flaring of the dipolar field lines, even in the absence of higher multipoles in the 
magnetic field (Scharlemann et al. 1978). This third effect is less important 14 than frame dragging if d{Bg/B, )d\nr ~ ovea /4 <C 
3k. This inequality translates into an upper bound on the surface polar magnetic field, 

/ P \ 3/7 

B NS < 1.3 x 10 16 ( — J G. (A32) 
The dipole fields of all magnetars easily satisfy this inequality. 



We consider the simplest case of an aligned rotator in making this estimate. 



